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The  active  intensity  and  wave  fronts  and  the  reactive 
intensity  measured  at  1000Hz  on  a  plane  perpendicular 
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measurements  were  made  in  reverberant  conditions . 
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The  potential  energy  and  resultant  phase  speed  measured 
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An  expanded  view  of  the  measurement  over  the  truncated 
sphere  top  plate  in  the  shadow  zone  of  a  rigid  disc. 

The  frequency  is  2570Hz  and  the  source  is  53cm  to  the 
left  of  the  disc  top . 
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top  plate . 


lax* 


%» , 


LIST  OF  TABLES 


Table 

6.1 

8.1 


The  real  and  imaginary  parts  of  the  impedance  normal 
to  the  baffle,  averaged  over  the  tube  opening.  The 
frequencies  of  200Hz  and  600Hz  ar  at  resonances  of  the  tube 
and  400Hz  is  at  the  tube  anti-resonance . 

A  summary  of  the  distinguishing  characteristics  of  the 
active  intensity,  reactive  intensity,  potential  energy, 
kinetic  energy,  and  the  real  and  imaginary  parts  of  the 
specific  acoustic  impedance  near  vibration  sources, 
resonators,  and  diffraction . 


<Oi  l? 


'.WW  wu 


BTuronufW 


WWW  UtlAJU-Jl  wfffwwwmraro  mu>  jiwihwwmilct 


A 


LIST  OF  SYMBOLS 


xviii 


amplitude  of  the  it;*  point  source 
speed  of  sound 
resultant  phase  speed 
separation  of  the  point  sources 
general  source  term  (complex  function) 

Heaviside  step  function 
active  intensity  vector 

complex  intensity  vector  (complex  function) 
instantaneous  intensity  vector 
instantaneous  intensity:  x-component 
instantaneous  intensity:  y-component 
unit  vector  in  the  x-direction 
unit  vector  in  the  y-direction 
square  root  of  —1 
free  field  wavenumber  (^) 

free  field  wave  vector  in  the  direction  of  wave  propagation 
resultant  wavenumber  (~) 

cp 

resultant  wave  vector  perpendicular  to  the  resultant  wave  fronts 

Lagrangian 

unit  vector 

acoustic  pressure 

acoustic  pressure  (complex  function) 
pressure  amplitude 
pressure  due  to  the  i**1  point  source 
reactive  intensity  vector 


I 


v\v 


N-Cv 


nnrjwv 


r 

r< 

U 

S 

SIL 

SPL 

SQL 

t 

T 

Tn 

u 

u 

Ua 

Uax 

UaV 

Hi 

Ur 

Urz 


U 


ry 

V 

Xi 

Vi 

2, 


XIX 

general  spatial  coordinate 

distance  from  the  \th  point  source  to  the  observation  point 

vector  from  the  ith  point  source  to  the  observation  point 

surface  area  enclosing  the  volume  V 

should  intensity  level  in  dB 

sound  pressure  level  in  dB 

reactive  intensity  level  in  dB 

time 

kinetic  energy  density  or  one  period  of  oscillation 
kinetic  energy  due  to  the  particle  velocity  in  the  direction  n 
particle  velocity 

particle  velocity  (complex  function) 

particle  velocity  in  phase  with  pressure 

particle  velocity  in  phase  with  pressure:  x-component 

particle  velocity  in  phase  with  pressure:  y-component 

particle  velocity  due  to  the  ith  point  source 

particle  velocity  in  phase  quadrature  with  pressure 

particle  velocity  in  phase  quadrature  with  pressure:  x-component 

particle  velocity  in  phase  quadrature  with  pressure:  y-component 

potential  energy  density  or  volume 

spatial  coordinate  of  the  ith  point  source 

spatial  coordinate  of  the  Ith  point  source 

spatial  coordinate  of  the  ith  point  source 

specific  acoustic  impedance  in  the  direction  h  (complex  function) 


*  m* 


Greek  and  other  symbols 


Q£t-  relative  phase  of  the  \th  point  source 

0  scalar  component  of  the  Helmholtz  decomposition  of  the  active 
intensity  vector 

Vnm  phase  term  for  the  n**  and  mth  point  sources 

0  angle  of  incidence  relative  to  the  norman  to  a  surface 
k  constant 

Hn  coordinate  (x,y,  or  z)  of  the  nth  point  source 

un  coordinate  (x,y,  or  z)  of  the  nth  point  source 

III  power  radiated  by  one  point  source 

II3  power  radiated  by  three  acoustically  coupled  point  sources 
p  air  mass  density 

<f>  pressure  phase:  a  function  of  space 

ip  vector  component  of  the  Helmholtz  decomposition  of  the  active 
intensity  vector 
u  radial  frequency 
(  )*  complex  conjugate 

(  )  signifies  a  complex  function 


sfOTwwwwri 


^  1 


ACKNOWLEDGEMENTS 


xxi 


I  sincerely  thank  the  following  people  for  their  guidance  and  friendship: 
Professor  Jiri  Tichy,  thesis  advisor,  and  Dr.  Courtney  Burroughs,  member 
of  the  committee,  who  provided  invaluable  advice,  insights,  and  support;  Dr. 
Sabih  Hayek  and  Dr.  Carter  L.  Ackerman,  members  of  the  committee,  for 
their  constructive  criticism;  Dr.  Anthony  Romano,  Dr.  Michel  Josserand, 
Dr.  Gary  Elko,  Karl  Grosh,  Jean-Francois  Degeorges,  and  Eugene  Taschuk 
who  were  always  willing  to  give  constructive  comments  and  encouragement;  and 
Barbara  Crocken,  administrative  aide  for  the  Graduate  Program  in  Acoustics, 
and  Carolyn  Smith,  department  secretary,  for  their  invaluable  assistance.  Above 
all  I  thank  my  parents,  Dr.  J.  Adin  Mann,  Jr.,  and  Arlene  E.  Mann  and  sisters, 
Elizabeth  and  Jennifer,  for  their  loving  guidance  and  support. 

This  work,  excluding  section  6.3,  was  supported  by  the  Naval  Sea  Systems 
Command. 

I  also  sincerely  thank  Dir.  u.  Prof.  Dr.  R.  Martin  and  Prof.  Dr.- 
Ing.  J.  Meyer  who  invited  me  to  be  a  guest  at  the  Physicalisch-Technische 
Bundesanstalt  (PTB)  in  Braunschweig,  Federal  Republic  of  Germany,  during 
the  fall  of  1987,  also  Dr.-Ing.  Hans-0  Finke  and  the  others  in  the  laboratory 
whose  skill,  patience,  and  great  hospitality  made  the  visit  invaluable.  The  work 
presented  in  section  6.3  was  supported  by  the  PTB. 


IT« 


INTRODUCTION  AND  OUTLINE 


1.1.  Introduction 

Energy  transfer  between  a  vibrating  body  and  the  surrounding  fluid  can  only 
be  throughly  understood  by  studying  the  nearfield  where  the  energy  transfer 
takes  place.  Likewise,  the  energy  propagation  where  diffraction  occurs  must  be 
examined  in  the  nearfield.  These  and  other  acoustic  phenomena  can  be  carefully 
investigated  using  acoustic  intensity  analysis.  Therefore,  the  principle  aspiration 
of  this  thesis  is  to  present  a  physical  interpretation  of  acoustic  intensity  and  to 
demonstrate  the  application  of  the  understanding  of  the  physical  interpretations 
of  acoustic  intensity  to  enhance  the  use  of  the  acoustic  intensity  technique  as  an 
analysis  tool  of  the  acoustic  nearfield. 

The  fundamental  issues  of  intensity  have  not  been  fully  developed  in  the 
open  literature  because  convenient  measurement  equipment  has  been  developed 
only  in  the  last  ten  years.  By  1975  the  two  microphone  technique  was  being 
used  with  analog  real-time  signal  processing.  The  inclusion  of  digital  signal 
processing  using  the  fast  Fourier  transform  (FFT),  first  shown  by  Fahy  (1977), 
greatly  aided  in  the  convenience  and  accuracy  of  acoustic  intensity  measurements. 
Elko  (19 84)  further  developed  the  technique  showing  that  the  complex  intensity 
vector,  energy  density,  and  instrument  and  field  dependent  bias  errors  are  easily 
calculated  from  the  FFT  of  the  pressure  signal  from  two  closely  spaced  pressure 
microphones. 

Acoustic  intensity  was  first  introduced  as  a  tool  to  measure  the  time  averaged 
power.  Limiting  acoustic  intensity  analysis  to  time  averaged  power  incorporates 
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only  a  small  fraction  of  the  available  information.  In  fact  contemporary 
interpretations  of  the  time  averaged  intensity  vector  axe  sometimes  in  error. 

Beyond  the  time  averaged  sound  power,  the  complex  acoustic  intensity  and 
related  analysis  provide  the  total  energy  propagation  and  wave  properties  in  a 
sound  field.  Many  of  these  issues  are  imbedded  in  articles  and  books  published 
before  1965.  However,  these  classical  concepts  of  sound  fields  have  not  been 
directly  linked  to  acoustic  intensity.  Researchers  are  only  now  beginning  to 
employ  fundamental  acoustic  intensity  concepts  necessary  to  understand  basic 
phenomena  that  occur  in  the  nearfield  of  vibrating  bodies  such  as  sources,  sinks, 
diffraction,  and  reflection. 

The  process  of  interpreting  the  active  intensity  vortex  played  an  important 
role  in  the  development  of  this  thesis  when  it  became  obvious  that  the  conven¬ 
tional  interpretation  of  the  active  intensity  lead  to  an  incorrect  understanding  of 
energy  flux  in  the  vortex.  Work  on  this  thesis  began  by  exploring  the  physical 
interpretation  of  the  intensity  vectors. 

1.2.  Goals  and  Outline  of  Thesis 

A  great  deal  of  initial  emphasis  is  placed  on  basic  interpretations  of  the 
acoustic  intensity  vectors.  The  ultimate  goals  of  this  work  are  to  demonstrate 
that  acoustic  intensity  is  a  powerful  analysis  tool  for  investigating  the  acoustic 
nearfield  of  sources  and  scatterers  and  interpreting  energy  transfer  between 
sources  and  pressure  fields  in  surrounding  fluids,  and  also  to  expand  the 
applications  of  the  intensity  analysis.  Based  on  fundamental  theories  of  sound 
fields,  the  acoustic  intensity  and  energy  density  and  flux  patterns  are  linked  to 
fundamental  concepts  of  acoustics.  The  measurements  which  are  presented  are 
chosen  to  demonstrate  basic  phenomena.  Only  aspects  of  each  measurement 
relevant  to  fundamental  concepts  of  acoustic  intensity  and  energy  flux  axe 
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discussed  so  that  the  logic  of  the  analysis  can  be  followed.  With  these  ends 
in  mind  the  thesis  is  organized  as  follows. 

Chapter  2  investigates  the  energy  flux  described  by  acoustic  intensity.  The 
time  dependent  intensity  is  compared  with  the  time  independent  active  and 
reactive  intensity.  The  comparison  focuses  on  the  description  of  time  dependent 
energy  flux  and  its  uniqueness. 

Chapter  3  links  fundamental  wave  properties  to  the  vector  characteristics  of 
•-he  active  and  reactive  intensity.  This  shows  that  the  active  intensity  represents 
the  propagating  wave  fronts  and  the  reactive  intensity  characterizes  the  pressure 
distribution.  The  phase  speed,  a  basic  description  of  a  wave,  is  calculated  from 
the  active  intensity  and  potential  energy. 

Chapter  4  contains  three  topics  which  are  understood  by  combining  the 
energy  and  wave  properties  discussed  in  Chapters  2  and  3.  These  topics  are  the 
group  velocity,  the  reactive  intensity,  and  the  active  intensity  vortex  formation. 

Chapter  5  briefly  describes  the  measurement  equipment  and  the  method  for 
measuring  acoustic  intensity  and  calculating  the  pressure  phase  from  the  time 
independent  active  intensity  and  potential  energy. 

Chapters  6  and  7  detail  measurements  which  illustrate  the  intensity  analysis 
methods  discussed  in  the  preceeding  chapters.  Simple  configurations  are 
investigated  which  model  some  basic  aspects  of  more  complicated  configurations. 
The  phenomena  which  are  focused  on  are  acoustic  sources,  absorbers,  diffraction, 
and  resonators. 

Chapter  6  concentrates  on  sources  and  passive  absorbers.  Measurements 
are  described  for  a  configuration  of  two  piston  sources  with  a  tube  between  the 
pistons.  Measurements  were  performed  at  the  one  quarter  wavelength  resonance, 
the  half  wavelength  anti-resonance,  and  three  quarter  wavelength  resonance  of 
the  tube.  The  second  example  compares  the  change  in  the  radiation  impedance 


and  output  power  for  a  standard  calibration  source  as  a  function  of  the  distance 
of  the  source  from  a  rigid  surface.  The  impedance  and  power  were  measured 
with  the  intensity  technique. 

Chapter  7  contains  three  examples  of  diffraction.  The  diffraction  at 
impedance  discontinuities  is  shown  around  a  loudspeaker  and  a  circular  disc. 
The  circular  disc  measurements  compare  the  diffraction  for  two  incident  wave 
configurations  and  are  made  with  a  rigid  and  an  absorptive  disc.  The  final 
example  is  measurements  made  above  a  truncated  sphere  in  the  shadow  zone  of 
a  rigid  disc.  The  analysis  includes  diagnostics  of  the  structural  vibrations  of  the 
truncated  sphere  induced  by  the  incident  sound. 

Chapter  8  contains  the  conclusions  and  recommendations  for  further  work. 
A  summary  of  the  analysis  capabilities  of  the  complete  intensity  measurements 
is  included  so  that  this  powerful  analysis  technique  can  be  applied  to  a  great 
variety  of  fundamental  and  applied  situations. 


5 


Chapter  2 


ACOUSTIC  INTENSITY  AS  A  DESCRIPTION 
OF  ENERGY  TRANSFER 


2.1.  Introduction 


The  fundamental  intent  of  acoustic  intensity  is  to  describe  the  transfer  of 
energy  in  an  acoustic  field.  To  this  end  extensive  literature  has  been  developed 
on  understanding  the  energy  transfer  process.  The  majority  of  work  done  with 
acoustic  intensity  focuses  on  the  time  averaged  intensity,  which  is  commonly 
called  the  active  intensity.  The  intent  of  this  chapter  is  to  derive  the  intensity 
quantities  in  terms  of  a  general  expression  for  the  acoustic  pressure.  The  focus 
will  be  placed  on  a  comparison  and  interpretation  of  the  energy  flux  as  described 
by  the  instantaneous  intensity,  active  intensity,  and  reactive  intensity. 

2.2.  Equations  for  Pressure,  Velocity,  and  Instantaneous  Intensity 

There  are  many  examples  in  the  literature  for  the  development  of  the 
formalism  for  acoustic  intensity  (Schultz,  Smith,  and  Malme  1964;  Pierce  1981; 
Elko  1984).  A  nice  introductory  summary  of  the  literature  is  given  by  Elko  (1984) 
who  develops  Equation  (2.1)  from  the  linear  acoustic  equations  in  the  form  of  a 
corollary  of  the  conservation  of  energy. 
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where  p  is  the  density,  p  is  the  acoustic  pressure,  u  is  the  magnitude  of  the 
particle  velocity,  u  is  the  particle  velocity  vector,  Ii(r,t)d=  p(r,t)  u(r,t )  is  the 
instantaneous  intensity,  and  n  is  a  unit  vector  normal  to  the  surface  5  enclosing 
the  volume  V.  This  equation  signifies  that  the  time  rate  of  change  of  the  total 


internal  energy  of  a  volume  is  equal  to  the  energy  flux  (instantaneous  intensity) 
through  the  volume  surface.  (Contributions  from  any  body  forces,  e.g.  gravity, 
and  viscous  dissipation  are  ignored.) 

Before  proceeding,  we  must  introduce  a  general  functional  form  to  be  used 
in  this  thesis  as  a  basis  to  study  a  general  three-dimensional  pressure  field.  The 
formula  is  given  and  discussed  in  Morse  and  Feshbach  (1953,  p.  145): 

p(r,t)  =  P(r)  fj<f>{r)  -  et]  ,  (2.2) 

where  r  is  a  general  spatial  variable  and  P(r)  and  <f>[r)  are  real.  (The  complex 
variables  are  underlined  throughout,  to  clearly  distinguish  them  from  the  real 
functions.)  This  form  will  be  discussed  in  greater  detail  in  Chapter  3.  For  a 
mono-chromatic  pressure  field  Equation  (2.2)  can  be  written  as 

p(r,t)  =  P(»y<w*-*r»  .  (2.3) 

In  this  form  the  spatial  amplitude  and  spatial  phase  distribution  are  explicitly 
described,  from  which  all  linear  acoustic  quantities  can  be  derived  (Pascal  and 
Lu  1984;  Mann,  Tichy,  and  Romano  1987). 

The  particle  velocity  is  calculated  using  Euler’s  equation  and  Equation  (2.3): 

u(r,t)  =  T-^-Vp(r,t) 

J(fP  (2-4) 

=  —  (P(r)V<£(r)  +  jVP(r))ej(wW(r))  . 

ujp 

Comparing  the  particle  velocity  with  pressure  explicitly  shows  that  there  are 
two  components  of  the  particle  velocity  (Kurtze  1968).  The  component  PV <j> 
is  in  phase  with  pressure  and  VP  is  ninety  degrees  out  of  phase  with  pressure. 
Consequently,  the  particle  velocity  can  be  written  in  a  form  to  emphasize  this 


u(r,t)  =  [U^r)  +  , 
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where 

Ua{r)  =  ~PV<t>  , 

7  (2-6) 

The  instantaneous  intensity,  /,-,  is  calculated  from  Equation  (2.1),  where  the 
form  of  the  particle  velocity  expressed  in  Equation  (2.4)  will  be  used.  One  must 
also  note  that  Equation  (2.1)  is  defined  for  real  quantities  thus  only  the  real  parts 
of  the  pressure  and  particle  velocity  in  equations  (2.3)  and  (2.4)  are  used: 

Ii{r,t)  =  Re{p{r,t)}  •  Pe{u(r,t)} 

=  ^P2(r)(V<£(r))cos2(wf  -  <f>)  (2.7) 

-  P[r)(V P(r))sin2(ut  -  <f>)  . 

Zup 

The  time  dependance  has  purposely  been  maintained  by  using  the  instanta¬ 
neous  intensity  so  that  the  energy  transfer  process  can  be  discussed  explicitly  in 
its  true  nature.  The  pressure  and  particle  velocity  are  phase  shifted  in  time,  but 
their  product  at  every  instant  is  the  energy  flux  from  the  point  in  space  where 
these  two  quantities  are  described. 

2.2.1.  Focus  on  One  Infinitesimal  Volume  of  Fluid 

A  short  detour  is  taken  from  the  formalism  to  focus  on  the  behavior  of  the 
fluid  described  in  the  previous  section.  We  consider  an  infinitesimal  mass  of  fluid 
as  it  undergoes  acoustic  motion. 

Morse  and  Ingard  (1968,  p.  227)  formulate  a  derivation  of  the  internal  energy 
in  an  acoustic  field  where  a  one-dimensional  duct  is  used  to  view  the  forces  and 
reactions  of  a  small  volume  of  the  fluid.  The  extension  to  three  dimensions  is 
very  useful.  Let  us  observe  an  element  of  fluid  with  three-dimensional  motion 
and  emphasize  cases  where  the  pressure  and  velocity  are  not  in  phase.  Since 
the  velocity  is  represented  by  two  vectors  in  space,  Ua  and  Ur,  there  is  a  plane 


which  contains  both  the  vectors.  Thus  the  general  three-dimensional  motion  at 
a  single  element  of  fluid  is  locally  two  dimensional. 

The  element  of  fluid  is  regarded  as  an  infinitesimal  mass  such  that  the 
particles  which  make  up  the  mass  move  together.  We  will  examine  the  pressure, 
volume,  velocity,  displacement,  and  acceleration  of  the  mass. 

The  discussion  will  not  be  confined  to  any  specific  source  configuration.  We 
are  only  interested  in  representing  a  situation  where  the  pressure  and  velocity 
are  phase  shifted  in  time.  The  following  equations  are  used: 

v  =  vo  —  vo cos(t)  , 

Jix(t)  =  eos(t)(UaXcos(t)  -  Urxsin[t))  ,  (2.8) 

Iiy{t)  =  cos  (f  )(I7ay  cos  (t)  —  Urvsin(t))  . 

where  v  is  the  instantaneous  volume  of  the  infinitesimal  mass,  vo  is  the  initial 
volume,  Ua  and  Ur  are  the  two  quadrature  components  of  the  particle  velocity, 
and  t  varies  from  0  to  2ir  to  represent  one  cycle  of  the  motion.  The  quantities 
Iix{t)  and  Iiy(t)  are  the  two  perpendicular  components  of  the  instantaneous 
intensity  vector. 

It  should  be  noted  that  the  motion  is  indeed  very  small,  a  result  emphasized 
by  the  acoustic  velocity  c  being  much  greater  than  the  particle  velocity.  Thus  the 
position  of  the  surrounding  masses  relative  to  the  mass  we  are  observing  remains 
essentially  constant  during  the  motion  of  the  mass. 

Time  histories  of  the  incremental  mass  for  two  different  values  of  Ua  and  Ur, 
are  plotted  in  Figures  2.1  and  2.2  for  one  complete  cycle  of  its  motion  around  the 
equilibrium  position.  The  initial  volume  is  represented  by  dashed  circles  and  the 
instantaneous  volume  is  represented  by  solid  circles.  The  path  is  represented  by 
a  thin  elliptical  line  with  the  instantaneous  volume  represented  at  one-twelfth 
period  increments  along  the  path.  At  each  position  the  instantaneous  intensity 
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Figure  2.1 


Figure  2.2 


The  path  of  an  infinitesimal  volume  during  one  period 
( Ua  —  i  ,  Ur  =  j).  The  dashed  circle  represents 
the  initial  volume,  the  solid  circle  represents  the 
instantaneous  volume,  and  the  arrow  represents  the 
instantaneous  intensity. 


1  _he  path  of  an  infinitesimal  volume  during  one  period 
[Ua  =  *(-5)  +  j(’25)  ,  Ur  =  i(.3)  +  ,;(.6)].  The  dashed 
circle  represents  the  initial  volume,  the  solid  circle 
represents  the  instantaneous  volume,  and  the  arrow 
represents  the  instantaneous  intensity. 
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direction  is  represented  by  an  arrow  whose  length  is  proportional  to  the  relative 
instantaneous  intensity  magnitude. 

In  Figure  2.1  Ua  and  UT  are  the  same  magnitude  but  perpendicular.  The 
resulting  particle  motion  is  a  circle.  In  Figure  2.2  Ua  =  t(0.5)  +  j(0.25)  and 
Ur  =  t(0.3)  +  j(0.6)  which  results  in  an  ellipse.  In  these  two  figures  there  are 
instants  labeled  A  to  F.  Each  letter  labels  two  instants  where  the  instantaneous 
intensity  points  in  a  given  direction  with  the  same  magnitude. 

Consider  the  times  labeled  by  A  in  Figure  2.1.  The  motion  of  the  mass 
at  the  top  position  is  to  the  left.  The  instantaneous  volume  is  greater  than 
the  equilibrium  volume  so  the  acoustic  pressure  is  negative.  The  instantaneous 
intensity,  the  product  of  the  acoustic  pressure  and  particle  velocity,  points  to  the 
right  which  at  this  instant  is  opposite  the  direction  of  motion.  When  the  mass  is 
at  the  bottom  position  the  velocity  is  to  the  right  and  the  volume  is  compressed 
so  that  the  acoustic  pressure  is  positive.  Thus  the  instantaneous  intensity  points 
to  the  right,  as  it  did  with  the  mass  in  the  top  position.  With  the  mass  at 
the  bottom  the  direction  of  the  instantaneous  intensity  and  the  particle  velocity 
are  the  same.  Despite  the  difference  in  the  direction  of  motion  for  these  two 
instances,  the  instantaneous  intensity  indicates  that  the  mass  is  doing  work  on 
the  same  neighboring  mass  to  the  right  at  both  times. 

As  seen  in  Equation  (2.7)  the  instantaneous  intensity  is  a  function  of  2 wt  thus 
oscillating  at  twice  the  mass  motion  frequency,  so  that  there  are  only  two  instants 
where  the  instantaneous  intensity  points  in  a  given  direction.  The  variation  of 
the  direction  of  the  instantaneous  intensity  vector  is  also  clearly  seen  in  Figures 
2.1  and  2.2.  This  observation  is  central  to  arguments  against  conventional  views 
that  the  energy  propagates  with  the  time  averaged  intensity,  /. 

Consider  a  case  where  Ua  and  Ur  are  in  the  same  direction  or  one  of  these 
is  zero,  the  path  of  the  mass  would  be  a  straight  line.  Consequently,  the 
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infinitesimal  mass  we  axe  viewing  would  only  transfer  energy  to  or  from  the  two 
adjacent  masses  along  the  motion  path.  This  corresponds  to  conventional  views. 
However,  counter  examples  are  easily  constructed.  For  if  Ua  and  UT  are  both 
non-zero  and  not  parallel,  there  are  only  two  instances  during  the  motion  cycle 
when  the  infinitesimal  mass  we  are  observing  transfers  energy  to  specific  adjacent 
masses.  In  other  words,  over  one  period,  the  mass  we  are  observing  is  having 
some  work  done  on  it  and  in  turn  is  doing  work  on  all  the  surrounding  masses, 
not  just  the  two  masses  along  the  direction  of  the  time  averaged  instantaneous 
intensity. 

The  view  that  the  time  averaged  intensity  sufficiently  describes  the  energy 
transfer  ignores  this  fact  and  says  that  the  mass  we  are  observing  only  transfers 
energy  to  one  neighboring  mass  which  is  oriented  in  the  direction  of  the 
time  averaged  intensity.  This  issue  can  be  more  thoroughly  addressed  when 
considering  a  large  region  of  an  acoustic  field  and  looking  at  some  specific  cases, 
which  will  be  done  in  section  2.5. 

2.3.  Time  Independent  Intensity  Vectors 

The  development  of  the  time  averaged  quantities  has  been  the  subject  of 
numerous  papers,  so  the  equations  are  only  briefly  developed  here.  In  this 
discussion  the  time  independent  intensity  vectors  will  be  derived  emphasizing 
their  relationship  to  the  general  form  of  the  acoustic  pressure,  which  is  expressed 
by  Equation  (2.2)  in  terms  of  the  pressure  amplitude  and  phase  spatial 
distributions.  These  intensity  vectors  will  be  linked  to  the  instantaneous  intensity 
and  the  energy  transfer  associated  with  the  intensity  vectors  will  be  discussed. 

The  interpretation  of  the  intensity  vectors  is  not  limited  to  energy  considera¬ 
tions,  and  often  the  wave  interpretations  are  more  significant  when  investigating 
an  acoustic  field.  The  wave  interpretations  will  be  discussed  in  detail  in  Chapter 
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3.  Here  the  discussion  will  be  confined  to  energy  interpretations  because  of 
the  distinction  which  must  be  made  between  the  time  dependent  instantaneous 
intensity  and  time  independent  active  and  reactive  intensity. 

2.3.1.  Development  of  Equations 

The  complex  intensity  vector  is  defined  as  (Pascal  1981;  Elko  1984): 

IM  ='  a  M  “  +j«M  ■  (2.9) 

The  spatial  components  of  the  pressure  and  particle  velocity  defined  in  Equations 


(2.3)  and  (2.4)  are 


P(r)  =  P(r)e-^{r)  , 

u(r)  =  — [P(r)V<^(r)  4- ;V?(r)]c-J>(r)  . 
up 


(2.10) 

(2.11) 


Substituting  Equations  (2.10)  and  (2.11)  into  Equation  (2.9)  results  in  the 


following  equation  for  the  complex  intensity: 


ZcM  =  —  . 


(2.12) 


The  real  and  imaginary  parts  of  Equation  (2.12)  define  the  two  vectors  which 
make  up  the  complex  intensity  vector.  The  vector  /,  defined  as  the  real  part  of 
Equation  (2.12),  will  be  called  the  active  intensity,  which  is 


The  imaginary  part  of  is  called  the  reactive  intensity: 


W  =  £pPW^r) 


(2.13) 


(2.14) 


The  origins  of  the  equation  for  the  complex  intensity  is  based  in  the  fact 
that  the  real  part  of  l_c  1S  the  time  average  of  the  instantaneous  intensity.  The 
instantaneous  intensity  is  defined  by  the  power  balance  equation  as  the  pressure 
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at  an  instant  of  time  times  the  particle  velocity  at  that  same  instant;  however,  the 
time  average  of  this  quantity  can  be  calculated  using  the  mathematical  complex 
forms  of  pressure  and  velocity.  This  mathematical  manipulation  also  allows  the 
reactive  intensity,  Q ,  to  be  calculated. 

2.3.2.  Energy  Transfer  Considerations 

Comparing  the  equations  for  the  active  and  reactive  intensity,  Equations 
(2.13)  and  (2.14),  and  the  instantaneous  intensity,  Equation  (2.7),  the  vectors 
I  and  Q  can  be  substituted  into  the  instantaneous  intensity.  Thus  the 
instantaneous  intensity  can  be  written  explicitly  in  terms  of  the  time  independent 
vectors  I  and  Q  as  follows: 

7t(r, t)  =  2l(r)cos2 (ut  —  <j>)  +  Q(r)sm2(wt  —  <j>)  .  (2.15) 

Expanding  cos2  (ut  —  (f>)  =  |  +  |cos2(wf  —  <j>)  gives  a  more  useful  form  of  Equation 

(2.15) : 

7,(r,  t)  =  7(r)  +  I{r)cos2{ijJt  —  <j>)  +  Q(r)sin2(ut  —  4>)  .  (2.16) 

This  form  of  the  instantaneous  intensity  clearly  shows  the  contribution  of  the 
active  and  reactive  intensity  to  the  energy  transfer  process.  The  instantaneous 
intensity  is  composed  of  the  term  7(r)  which  is  independent  of  time.  The 
two  additional  terms  containing  the  active  and  reactive  intensity  oscillate  in 
quadrature  at  twice  the  frequency,  as  illustrated  in  Section  2.2.1. 

Conventionally,  the  instantaneous  intensity  is  time  averaged,  which  Equation 

(2.16)  clearly  shows  leaves  the  active  intensity.  However,  the  instantaneous 
propagation  of  energy,  described  by  the  instantaneous  intensity,  also  includes 
the  reactive  intensity.  This  is  evident  in  Figures  2.1  and  2.2  where  at  only  two 
instances  in  one  cycle  does  the  instantaneous  intensity  point  in  a  given  direction. 
These  two  figures  show  that  the  direction  of  the  energy  flux  is  biased  in  one 
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direction  by  the  active  intensity  term,  but  in  general,  the  energy  flux  oscillates 
around  the  direction  of  the  active  intensity.  The  active  intensity  does  describe  the 
average  direction  and  magnitude  of  the  instantaneous  intensity.  But  due  to  the 
averaging  process,  the  active  intensity  does  not  describe  the  total  energy  transfer, 
because  of  the  two  terms  in  Equation  (2.16)  which  are  ignored.  Quite  simply, 
the  average  direction  of  a  vector  does  not  necessarily  indicate  the  direction  of  the 
vector  at  all  times.  The  active  intensity  in  general  only  describes  the  energy  flux 
at  two  instances  in  one  cycle,  not  the  time  dependent  process  of  energy  transfer. 

Justifications  that  the  time  independent  active  intensity  adequately  describes 
energy  flow  often  involve  the  rationalization  that  the  reactive  intensity  although 
present,  is  of  no  consequence  because  its  contribution  time  averages  to  zero. 
Equation  (2.16)  contains  a  term  of  the  active  intensity  which  likewise  time 
averages  to  zero.  Consequently,  if  the  active  intensity  alone  is  used  to  describe 
the  energy  flux  over  all  time  then  there  is  a  contribution  of  the  active  intensity, 
I[r)cos2{<jJt  —  d>(r))  which  is  likewise  being  ignored  along  with  the  reactive 
intensity. 

Interpretations  of  Equation  (2.16)  is  the  basis  of  the  argument  against  the 
conventional  view  of  active  intensity  as  a  total  description  of  energy  transfer. 
The  active  intensity  has  great  significance,  but  I  only  describes  the  net  flux  of 
energy  at  a  point,  I  does  not  indicate  in  general  where  the  energy  goes  or  how 
it  arrives  at  the  point. 

Take  the  example  of  a  pure  standing  wave,  where  the  active  intensity  is 
zero.  Conventional  viewpoints  would  indicate  that  there  is  no  energy  flux  in 
the  standing  wave.  However,  the  reactive  intensity  is  non-zero.  Therefore,  from 
Equation  (2.16),  the  instantaneous  intensity  is  non-zero,  indicating  that  there 
is  energy  propagating  in  the  standing  wave.  At  every  point  the  energy  flux 
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will  time  average  to  zero,  but  the  existence  of  the  acoustic  field  requires  energy 
propagation. 

These  arguments  go  beyond  semantics.  It  is  important  that  the  proper 
physical  meanings  are  given  to  the  instantaneous  intensity,  active  intensity,  and 
reactive  intensity  vectors,  in  terms  of  energy  propagation.  These  views  do  not 
diminish  the  significance  of  the  vectors,  but  try  to  specify  the  insight  the  vectors 
provide  for  the  energy  propagation  in  an  acoustic  field. 

2.4.  Sources  of  Acoustic  Energy 

The  net  outflow  of  a  vector  is  the  divergence  of  the  vector.  Since  derivations 
of  the  divergence  of  I  and  Q  are  well  documented,  (Schultz,  Smith  and  Malme 
1964;  Elko  1984)  only  the  results  will  be  stated  here: 


V  •  /  =  0  , 


V  •  Q  =  -2w(T  ~V)  =  -2 wL  , 


(2.17) 


(2.18) 


where  L,  the  Lagrangian,  is  the  difference  of  the  potential  and  kinetic  energy. 

2.4.1.  Sources  of  Active  Intensity 

Equation  (2.17)  is  of  great  interest  because  it  states  that  the  divergence  of  the 
active  intensity  is  zero  everywhere  it  is  defined.  However,  this  term  is  contained 
in  the  conservation  of  energy  corollary  stated  in  Equation  (2.1),  which  in  time 
averaged  form  is 

/ JIM + ^iv  -  -  III v  • Idv  •  (2-i9) 

This  equation  is  the  basis  of  power  measurements  which  relate  the  power  radiated 
from  a  source  in  volume  V  to  the  active  intensity.  However,  substituting  Equation 
(2.17)  into  Equation  (2.19)  the  time  averaged  change  in  energy  in  the  volume  V 
is  zero.  Consequently,  if  the  source  which  causes  the  acoustic  field  in  volume  V 
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does  radiate  power,  the  source  is  not  in  the  volume  V .  Therefore  we  must  be 
careful  in  constructing  the  volume  in  which  the  intensity  is  defined. 

The  conservation  of  energy  corollary,  Equation  (2.1),  is  derived  for  a  source 
free  region.  Thus  the  intensity  vector,  for  our  use,  is  only  defined  in  regions  of 
no  sources.  If  the  acoustic  field  around  a  source  is  of  interest  then  the  volume 
sketched  in  Figure  2.3  can  be  constructed.  The  volume  integral  can  subsequently 
be  converted  to  a  surface  integral  using  Green’s  Theorem.  However,  the  surface 
integral  is  not  only  over  54,  which  is  conventionally  used,  but  also  includes  Si, 
S2,  and  S3.  These  added  surfaces  allow  the  construction  of  the  volume  around 
but  not  including  the  sources.  The  surface  integral  will  be 

V  •  IdV 


ill  V  '  UV  -  fls,  1  ■  +  J is.  7 '  iSl  +  I L  7 '  ^ 


//, 


(2.20) 


I -dS4  . 


The  paths  S2  and  S3  will  give  the  same  result  but  opposite  in  sign,  therefore 
cancel.  The  integral  over  Si  and  54  are  left  and  must  cancel,  which  means  that 
the  energy  entering  V  through  Si  leaves  through  S4  ,  and  that  the  power  output 
by  the  source  is  quantified  by  the  integral  over  Sj  or  S4. 

This  observation  produces  the  same  result  which  has  always  been  used  for 
power  measurements,  but  it  is  important  to  realize  the  formality  of  this  argument 
which  is  based  on  the  fact  that  the  intensity  vectors  are  only  defined  in  source 
free  regions. 

Pascal  (1985)  derived  a  Helmholtz  decomposition  of  the  vector  /  which  split 
it  into  a  part  defined  from  a  gradient  of  a  scalar  and  another  part  from  the  curl 
of  a  vector: 

/  =  Vi3  +  Vxt?  .  (2.21) 

The  divergence  of  the  term  V  x  xjj  is  zero  by  definition  of  the  divergence  and 
curl  operators.  However,  the  divergence  of  the  term  V/?  is  in  general  not 
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zero.  Consequently,  Pascal  concluded  from  Equation  (2.19)  that  V/3  accounts 
for  radiated  power.  But,  as  already  stated,  V  •  /  =  0,  thus  V  •  V/3  must  be  zero. 
In  light  of  the  previous  discussion  of  the  volume  in  which  intensity  is  considered, 
the  conclusion  that  V/3  does  not  account  for  the  time  averaged  output  power  of 
a  source  is  valid. 

This  same  conclusion  can  be  reached  by  considering  the  case  of  the  radiation 
of  N  point  sources,  where  the  total  pressure  and  particle  velocity  are 


P(r)  =  J2  E,(r)  > 


5(r)  =  ^u,(r)  , 


where 


(r)  = 

c‘  '  rt- 
u-M  = 

jup  rj1 


(2.22) 


(2.23) 


(2.24) 


(2.25) 


The  constants  A,  and  a,-  are  the  relative  amplitude  and  phase  of  the  i-th  source, 
and  r,  and  f*  are  the  distance  and  direction  vector  from  the  i-th  source  to  the 
measurement  point  r: 


r 2  =  (x  -  Xi)2  +  (y  -  y,)2  4-  {x  -  a,)2  , 


rx  =  (x  -  Xi,  y  -  y,-,  a  -  2.)  . 


(2.26) 


(2.27) 


The  complex  intensity  vector  7e  can  be  calculated  from  equations  (2.22)  and 
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(2.23): 


lc(r)  =  \m 


=  fL\±AiL^l] 

2jujpI^  n  JL  fr{  rf  J 


*  2  upr*  1 


(2.28) 


TV  TV 


+  jjr  y*  iblAjf  e-m^-rm)-(an-a, 

2a/p  rnr2,  m 


2a/a  rnr3 

n-l  m  =  !  '"'m 

m  ?6  n 

Since  / c  =  /  +  jQ  and  defining  r?nm  =  fc(rn  -  rm)  —  (a„  —  am)  Equation  (2.28) 
can  be  separated,  such  that 
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_  1  T‘  1  _  ^  f 

I{r)  —  o j  -  ^  >  Ai  ^3  ~h  ^  }  !  )*  ^  ^mAn  ~2~{krmCOsr)nm  —  sinr)nm}  ,  (2.29) 
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=  2wp  ^-'j4,2H'+2wo  51  51  A»xAn^-^-{cosrfnm-krmsmr/nm)  .  (2. 

t=l  »  K  n=  1  m=  i  n  rn 

m  ^  n 

Now  calculate  the  curl  of  the  intensity  defined  in  Equation  (2.29): 


TV  TV 


Vx7(r)  =  ^S  ^  ^{[(fcVm  +  l)«n>?f 


n= 1 m=n  + 1 


+  k{rm  ~  rn)cosr,nm  <  tynzm,  eznzm,  ezym  > 


,  (2.31) 


where  <  £j(rr,»  enm»  £«m  >  's  a  vector  which  has  the  components  =  ix{v>m  — 
yn)  +  T/(pn  —  pm)  ;  p  =  i,  y,  or  2,  u  =  x,  y,  or  z.  Equation  (2.31)  shows  that  the 
part  of  the  equation  for  the  active  intensity  with  the  cross  terms  proportional  to 
AmAn(n  7^  m)  is  the  part  of  I  which  has  a  non-zero  curl.  Consequently,  for  a 
distribution  of  point  sources  the  active  intensity  can  be  split  as  Equation  (2.21) 


in  the  following  manner: 


V0 -*][>?£  , 


(2.32) 


Eli 


KW 


20 


JS.I‘LfUVI'1 


AT  JV 


xjj  =  k  EE  ■^m^n  2  (.^TmCOSTlnm  "T"  SlTlTlnrn ) 

n  =  l  -"  =  1  rnrm 


(2.33) 


Equation  (2.32)  means  that  V/3  is  the  vector  sum  of  the  active  intensity  from 
each  point  source  ignoring  the  coupling  of  the  sources.  Further,  V  x  rp  defined 
in  Equation  (2.33)  represents  the  effect  of  coupling  between  the  point  sources. 

As  an  example  we  will  calculate  the  total  power  II3  radiated  by  three  point 
sources  compared  to  the  power  TTi  radiated  by  one  point  source.  The  three  point 
sources  are  in  a  line  separated  by  distance  d  with  relative  amplitudes  +1,  —1,  +1. 
Using  the  expressions  in  Skudrzyk  (1971), 


II3  sinkd  sin2kd 

—  =  3  -  4 - 1 - 

nx  kd  kd 


(2.34) 


Figure  2.4  represents  Equation  (2.34)  as  a  function  of  kd,  clearly  shows  the  effect 
of  mutual  coupling  of  the  sources  especially  at  low  kd.  To  identify  the  terms  in 
Equation  (2.34)  the  power  can  be  calculated  from  the  intensity  by  integrating 
I  over  a  surface  enclosing  the  sources.  Without  performing  the  integration 
explicitly,  the  power  associated  with  V/?  is^clearly  3IT 1 .  Consequently,  the  power 
associated  with  V  x  xjj  is  —4  -^d  +  Therefore,  in  the  region  excluding  the 

sources,  V/3  does  not  totally  account  for  the  energy  propagating  to  the  farfield, 
because  the  mentioned  terms  due  to  the  source  coupling  are  not  negligible. 

2.4.2.  Sources  of  Reactive  Intensity 

A  very  nice  development  of  the  sources  of  reactive  intensity  is  presented  by 
Elko  (1985).  He  derives  an  expression  from  the  complex  power  and  the  divergence 
of  the  complex  intensity  vector,  which  is  much  the  same  as  the  power  balance 


for  the  active  intensity: 


JJ  Q-ndS  =  11/ +  ///  2tt LdV  , 


(2.35) 
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Figure  2.4 


The  ratio  of  the  power  radiated  by  three  sources  to 
the  power  radiated  by  one  source.  The  sources  are 
separated  by  18cm  with  the  center  source  out  of  phase. 
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where  II;  is  the  imaginary  part  of  the  complex  power  radiated  by  a  source,  L  is 
the  Lagrangian,  which  enters  from  the  divergence  of  the  reactive  intensity,  V  is 
the  volume  enclosing  a  source,  and  S  is  the  surface  area  of  the  volume. 

Equation  (2.35)  shows  that  the  reactive  energy  enters  the  volume  through 
its  boundaries  (sources)  as  signified  by  the  term  IT/,  but  there  is  also  the 
energy  which  is  originating  in  the  volume,  represented  by  the  term  including  the 
Lagrangian.  The  integral  of  the  Lagrangian  signifies  that  the  reactive  energy  is 
present  because  of  a  difference  of  the  kinetic  and  potential  energy.  It  is  interesting 
that  it  includes  the  difference  of  the  total  kinetic  energy  which  has  a  component 
associated  with  the  particle  velocity  accounting  for  the  active  intensity. 

The  reactive  intensity  exists  if  there  is  a  nonvanishing  pressure  gradient. 
With  a  nonvanishing  pressure  gradient  there  is  an  imbalance  in  the  potential 
and  kinetic  energies.  Thus,  there  is  energy  which  is  removed  from  the  potential 
energy  density  which  is  needed  for  maintaining  the  pressure  gradient.  This  is 
energy  which  is  input  to  a  region  during  one-half  of  a  cycle  then  work  exerted 
back  by  the  region  during  the  other  one-half  of  a  cycle.  Consequently,  energy  is 
removed  from  that  energy  which  is  radiated  to  infinity  by  a  source  to  maintain 
a  nonvanishing  pressure  gradient  associated  with  the  imbalance  in  the  potential 
and  kinetic  energies. 

2.5.  Instantaneous  Energy  Flux  Lines 

As  discussed  in  Section  2.3.2.,  the  instantaneous  intensity  describes  energy 
propagation  through  a  point.  To  follow  energy  through  space,  lines  must  be 
drawn  parallel  to  the  instantaneous  intensity  vector.  The  resulting  instantaneous 
intensity  flux  lines  will  create  a  complicated  picture  because  they  are  a  function 
of  space  and  time.  These  lines  clearly  show  the  distinction  between  the  time 


dependent  description  of  energy  transfer  (instantaneous  intensity)  and  the  time 
independent  description  of  power  flux  (active  intensity). 

The  scheme  to  construct  the  instantaneous  intensity  flux  lines  is  represented 
in  Figure  2.5  for  a  hypothetical  example.  The  instantaneous  intensity  vectors 
at  six  points  are  drawn  for  two  times,  to  and  to  +  At.  The  resulting  flux  line  is 
drawn  with  a  dashed  line.  The  variation  of  the  instantaneous  intensity  vectors 
with  the  increment  of  time  is  greatly  exaggerated  in  Figure  2.5,  as  is  the  distance 
between  the  points  between  which  the  flux  line  is  drawn. 

The  instantaneous  intensity  flux  line  begins  in  Figure  2.5a  at  point  1  at  time 
t0.  The  flux  line  is  drawn  in  the  direction  of  the  instantaneous  intensity  at  time 
*o  at  point  1,  which  is  toward  point  5.  Thus  the  energy  propagates  through 
point  1  toward  point  5  at  speed  c.  The  propagation  distance  is  A r,  so  that  this 
propagation  takes  time  At  =  Ar.  At  time  to  +  At  the  instantaneous  intensity 
vectors  at  the  six  points  are  different,  as  represented  in  Figure  2.5b.  The  energy 
now  propagates  through  point  5  toward  point  3  in  time  At. 

The  instantaneous  intensity  flux  line  for  this  example  is  drawn  from  point  1 
to  point  5  to  point  3  and  continues  through  space  as  time  increases.  The  resulting 
line  represents  the  energy  propagation  as  a  function  of  space  and  time,  although 
only  the  spatial  dependence  can  be  shown.  The  direction  of  the  instantaneous 
intensity  at  point  1,  the  initial  point,  varies  over  one  half  of  a  cycle.  When  the 
line  begins  at  some  other  time  at  point  1  the  path  will  begin  in  another  direction 
and  continue  on  another  path.  Only  at  times  separated  by  one  half  of  a  period 
does  the  line  beginning  at  point  1  trace  the  same  path.  Consequently,  there  is 
no  one  line  which  describes  the  energy  propagation  at  all  times. 
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a)  time  to 


b)  time  t0  +  A t 


Figure  2.5 


Schematic  of  constructing  instantaneous  intensity  flux 
lines.  The  instantaneous  intensity  vectors  are  repre¬ 
sented  at  six  points  at  times  to  and  to  +  At,  and  the 
resulting  flux  line  (dashed)  is  drawn. 
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2.5.1.  Instantaneous  Intensity  Flux  Lines  for  Three  Point  Sources 

Consider  the  field  radiated  by  point  sources  arranged  on  a  line  with  18cm 
center  to  center  spacing.  The  two  outer  sources  are  in  phase  and  the  center  source 
vibrates  in  opposite  phase.  This  source  configuration  is  used  because  it  provides 
examples  of  high  mutual  coupling  and  low  mutual  coupling.  Two  frequencies 
will  consequently  be  discussed:  (l)  2070  Hz  (low  mutual  coupling,  the  sources 
are  separated  by  1.1  A)  and  (2)  200  Hz  (high  mutual  coupling,  the  sources  are 
separated  by  0.1A).  The  active  and  reactive  intensity  for  these  frequencies  are 
plotted  in  Figures  2.6  and  2.7.  These  are  commonly  discussed  examples  so  the 
analysis  will  be  limited. 

At  2070  Hz,  with  the  low  mutual  coupling,  the  active  and  reactive  intensity 
point  out  of  the  three  sources.  Thus  one  can  conclude  that  all  the  energy  in  the 
region  propagates  away  from  the  sources.  At  200  Hz,  high  mutual  coupling,  the 
active  intensity  points  into  the  center  source,  leading  to  the  conclusion  that  the 
center  source  absorbs  power.  However,  the  reactive  intensity  points  out  of  each 
source  since  the  pressure  is  a  local  maximum  above  the  sources,  independent  of 
the  direction  of  the  active  intensity. 

Now  look  at  the  instantaneous  intensity  flux  lines;  Figure  2.8  and  2.9.  The 
energy  is  followed  from  four  initial  points  marked  by  circles.  Two  initial  times  are 
chosen  so  there  are  two  paths  which  leave  each  point.  The  direction  of  the  energy 
propagation  is  indicated  by  an  arrow.  There  is  an  infinite  set  of  paths  leaving 
each  point  at  different  instances  during  one  period;  however,  for  the  purposes 
here  two  initial  times  adequately  demonstrate  an  important  point. 

Figure  2.8  is  intuitively  very  easy  to  understand.  Each  path  shows 
propagation  away  from  the  three  sources.  This  would  be  expected  because  the 
sources  have  low  coupling  and  the  active  and  reactive  intensity  vectors  point 
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Figure  2.6  The  calculated  active  intensity  and  reactive  intensity 
around  three  point  sources  separated  by  18cm  with  the 
center  source  out  of  phase.  The  frequency  is  2070Hz. 
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Figure  2.7 


The  calculated  active  intensity  and  reactive  intensity 
around  three  point  sources  separated  by  18cm  with  the 
center  source  out  of  phase.  The  frequency  is  200Hz. 
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Figure  2.8  Instantaneous  intensity  flux  lines  at  2070Hz.  (Three 
point  sources  with  the  center  source  out  of  phase.) 
The  circles  represent  the  starting  positions  of  the  lines. 
Two  lines  leave  each  starting  position  one  at  times  tQ 
and  t0  + 
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Figure  2.9  Instantaneous  intensity  flux  lines  at  200Hz.  (Three 
point  sources  with  the  center  source  out  of  phase.) 
The  circles  represent  the  starting  positions  of  the  lines. 
Two  lines  leave  each  starting  position  one  at  times  t0 
and  ta  -f  7T. 
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out  of  the  sources.  Thus  the  present  views  that  the  active  intensity  show 
energy  propagation  are  essentially  true  since  it  would  conclude  that  the  energy 
propagates  from  the  sources. 

However,  Figure  2.9  shows  a  rather  surprising  result,  because  each  path, 
although  erratic,  does  go  to  the  farfield.  Thus  energy  at  each  of  the  four  starting 
positions  does  get  to  the  farfield.  This  directly  contradicts  present  views  that  the 
active  intensity  shows  the  path  of  energy  propagation,  since  for  this  frequency, 
Figure  2.7,  the  active  intensity  would  say  that  all  energy  in  the  four  starting 
positions,  shown  in  Figure  2.9,  goes  into  the  center  source.  Consequently  it  is 
invalid  to  say  that  the  active  intensity  shows  the  path  of  energy  propagation, 
because  it  is  a  time  averaged  quantity  which  does  not  include  the  contribution 
of  the  reactive  intensity  to  the  energy  transfer. 

It  was  stated  that  in  each  period  of  the  acoustic  motion  an  infinite  set 
of  different  paths  from  each  point  can  be  constructed.  Consequently  an 
important  issue  is  the  practicality  of.studying  simultaneously  the  time  and  spatial 
dependance  energy  propagation.  It  is  obvious  that  a  proper  description  of  energy 
propagation  must  include  the  reactive  intensity  along  with  the  active  intensity; 
however,  this  is  an  impractical  view  because  there  is  not  one  path  along  which 
the  energy  propagates  at  all  time.  Thus  one  must  turn  to  the  time  independent 
vectors,  I  and  Q,  and  remember  that  the  active  intensity  vector  quantifies  in 
magnitude  and  direction  the  time  averaged  power  flux,  not  the  path  of  energy 
propagation.  It  should  also  be  understood  that  the  energy  in  the  reactive 
intensity  does  not  simply  oscillate  back  and  forth  but  is  indeed  transferred  to 
the  farfield  in  the  form  of  the  energy  propagation  described  by  the  instantaneous 
intensity. 


A  few  authors  have  published  papers  about  constructing  lines  which  they  call 
energy  streamlines  that  are  drawn  tangent  to  the  active  intensity  vector  (Wa¬ 
terhouse,  Yates,  Feit,  and  Liu  1985;  Skelton  and  Waterhouse  1986;  Waterhouse 
and  Feit  1986;  Waterhouse  1987).  The  claim  is  that  these  lines  show  the  flow  of 
energy  through  space.  However,  as  already  discussed,  the  instantaneous  intensity 
is  the  description  of  energy  flow.  Consequently,  to  follow  energy  through  space, 
energy  flux  lines  must  be  drawn  parallel  to  the  instantaneous  intensity.  The 
resulting  instantaneous  intensity  flux  lines  show  the  distinction  between  the  time 
dependent  process  of  energy  transfer  and  the  time  averaged  power  flux. 


Chapter  3 


WAVE  PROPERTIES  DESCRIBED  BY  INTENSITY 

3.1.  Introduction 

Intensity  is  used  primarily  to  describe  the  transfer  of  energy  in  an  acoustic 
field.  However,  the  vector  properties  of  the  active  and  reactive  intensity  are 
linked  directly  to  properties  of  wave  forms.  These  interpretations  begin  from  the 
general  description  of  acoustic  pressure  which  was  used  to  develop  the  energy 
transfer  characteristics  in  the  previous  chapter.  The  result  of  this  analysis  is 
a  method  of  representing  the  wave  structure.  This  proves  especially  useful  in 
the  nearfield  where  the  acoustic  field  is  normally  viewed  as  very  complicated. 
The  visualization  of  the  nearfield  which  is  provided  by  the  active  and  reactive 
intensity  is  a  valuable  use  of  the  acoustic  intensity  analysis  along  with  the  energy 
interpretations. 

There  have  been  many  authors  who  have  addressed  the  issues  of  wave 
propagation.  Commonly  the  wave  structure  is  discussed  in  terms  of  the  energy 
density,  wave  shape,  phase  speed  and  group  velocity.  The  new  contribution  of 
this  chapter  is  to  link  these  well  established  ideas  to  the  intensity  analysis. 

3.2.  General  Functional  Form  of  Acoustic  Pressure 

The  essence  of  this  chapter  is  the  interpretation  of  the  general  functional  form 
for  the  acoustic  pressure  which  was  briefly  introduced  in  the  previous  chapter. 
The  pressure  is  described  in  three  dimensions  by  the  following  equation: 

p{r,t)  =  P(r)  f[<f>(r)  -  ct]  ,  (3.1) 

where  r  is  a  general  coordinate,  P(r )  is  a  positive  real  function  of  the  general 
coordinate,  and  4>(r )  is  a  real  function  of  the  general  coordinate.  This  form 
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clearly  separates  the  pressure  into  two  parts:  (1)  the  pressure  amplitude  spatial 
distribution  specified  by  P(r )  and  (2)  the  phase  spatial  distribution  specified  by 
-  ct]. 

An  important  aspect  of  this  description  is  the  wave  shape  and  wave 
propagation  described  by  f[4>{r )  —  ct\.  Contained  in  this  function  is  the  wave 
front  which  is  inherent  in  our  concepts  of  acoustic  fields.  For  this  thesis  the 
function  will  be  studied  for  monochromatic  fields,  so  that  the  pressure  is  written 
as 

p(r,t)  =  .  (3.2) 

In  classical  theories  waves  are  defined  as  surfaces  of  constant  phase,  thus  this 
functional  form,  although  specifying  the  time  dependence  as  harmonic,  maintains 
a  general  wave  front  shape  defined  by  surfaces  where  d>(r)  =  constant. 

This  functional  form  is  not  the  most  general  representation  of  wave  propaga¬ 
tion.  Whitham  (i960;  1961)  discusses  the  notion  of  a  wave  form  from  its  transient 
development  to  the  steady  state;  consequently,  addressing  the  time  and  spatial 
variation  of  an  acoustic  wave  form.  However,  for  the  purpose  of  this  thesis  it 
is  sufficient  to  investigate  a  steady  state  field.  These  fields  can  be  decomposed 
into  frequency  components,  so  it  is  also  valid  to  discuss  monochromatic  fields. 
Consequently  the  developments  can  be  focused  on  the  consequence  of  a  general 
wave  shape  4>{r )  with  a  pressure  distribution  P[r )  recognizing  this  as  a  step  to 
extending  the  understanding  toward  broad  band  and  transient  phenomena. 

3.3.  Differential  Equations  for  P{r )  and  <fr{r) 

A  constraint  on  the  generalized  form  of  the  pressure  is  that  it  satisfies  the 
wave  equation.  The  function  P(r )  and  cf>(r)  can  not  be  chosen  arbitrarily.  But 
differential  equations  can  be  developed  which  relate  the  two  functions  so  that 
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p(r,  t)  satisfies  the  wave  equation.  This  interdependence  is  derived  by  substituting 
the  pressure  into  the  Helmholtz  equation: 


v2p(r>*)  -^2~  =£(r,0  -  (3.3) 

where  g{r,t)  is  a  source  term.  However,  the  intensity  is  developed  for  a  source 
free  region  so  the  source  term  g(r,  t)  is  zero.  Substituting  Equation  (3.2)  into 
Equation  (3.3)  gives 

[V2P-PV<£-V<£-j(2VP.V<£  +  pV2^)  +  fc2P]e^wt-^  =  0  .  (3.4) 

The  real  and  imaginary  parts  of  this  equation  can  be  separated  so  that  the 
following  differential  equations  can  be  written: 

V2P  +  [k2  -  V<f>  ■  V<£]P  =  0  ,  (3.5) 

PV  V  +  2V<£  •  VP  =  0  ,  (3.6) 

where  k  =  which  is  independent  of  position.  Equations  (3.5)  and  (3  6) 
are  non-linear  partial  differential  which  couple  the  phase  characteristics  of  the 
pressure  to  the  pressure  amplitude  variation. 

3.4.  Vector  Characteristics  of  /  and  Q 

Before  going  any  further  we  need  to  attach  a  physical  meaning  to  the  vector 
characteristics  of  the  active  and  reactive  intensity.  From  the  functional  form  of 
the  acoustic  pressure,  Equation  (3.2),  the  following  equations  were  developed  in 
the  previous  chapter  for  I  and  Q: 

fy)  i^pJ(r)v*(r)  ■  (3-7) 

CM  5^p(r)vp(r)  '  (3-8) 
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The  active  intensity  has  the  direction  of  the  spatial  gradient  of  the  pressure 
phase,  consequently  /  is  perpendicular  to  surfaces  of  constant  phase.  By 
definition,  the  surfaces  of  constant  phase  are  wave  fronts;  so  the  active  intensity 
is  perpendicular  to  the  resultant  wave  fronts  in  an  acoustic  field.  The  term 
“resultant”  is  used  because  these  wave  fronts  are  a  result  of  the  interference  of 
many  waves.  For  instance,  with  a  plane  wave  incident  on  a  reflecting  boundary, 
the  acoustic  field  can  be  considered  as  the  interference  of  the  direct  plane  wave 
and  the  reflected  plane  wave.  This  interference  produces  the  resultant  waves  to 
which  the  active  intensity  is  perpendicular. 

As  noted  by  numerous  authors  the  reactive  intensity  is  proportional  to  the 
negative  of  the  spatial  gradient  of  the  pressure  amplitude  (Elko  1984;  Tichy 
1984).  Consequently,  Q  points  normal  to  surfaces  of  constant  pressure.  More 
specifically,  Q  points  out  of  pressure  maxima  and  into  pressure  minima.  This 
proves  very  useful  for  visualizing  the  fine  structure  of  the  pressure  distribution. 

With  the  interpretations  of  the  active  and  reactive  intensity  the  acoustic  field 
can  be  organized  into  resultant  wave  fronts  and  a  pressure  distribution.  The 
idea  of  organizing  the  acoustic  field  is  central  to  the  powerful  insight  the  acoustic 
intensity  analysis  provides  for  visualizing  and  interpreting  the  nearfield  of  sources 
and  scatterers.  This  will  become  evident  with  several  examples  in  Chapters  6 
and  7  where  these  two  simple  ideas  are  applied  to  otherwise  complicated  acoustic 
nearfields. 

3.5.  Generalized  Phase  Speed 

The  active  intensity  has  been  shown  to  visualize  the  resultant  wave  fronts. 
Implicit  to  the  concept  of  a  wave  front  is  the  speed  of  propagation.  Convention¬ 
ally,  the  speed  of  sound  c,  is  associated  with  the  propagation  speed  of  simple 
waves  such  as  plane  waves  and  spherical  waves.  However,  in  general  the  speed  of 


sound,  c,  is  different  than  the  speed  of  wave  front  propagation,  cp.  For  instance, 
this  is  recognized  in  a  waveguide  where  the  phase  speed  is  not  c.  The  extension 
to  the  general  acoustic  field  has  been  done  in  the  historic  literature  (Rayleigh 
1877,  Lamb  1904),  but  is  included  here  to  link  the  phase  speed  with  the  active 
intensity. 

The  general  phase  speed  for  the  surfaces  of  constant  phase  is  calculated  by 
setting  the  pressure  phase  equal  to  a  constant  for  a  change  in  time  and  position 
normal  to  the  surfaces  of  constant  phase.  Thus, 

(ut0  -  <f>(r0))  =  [u(t0  +  At)  -  <f>(r0  +  Ar)]  .  (3.9) 

This  is  equivalent  to  setting  the  total  derivative  of  the  pressure  phase  with  respect 
to  time  and  distance  normal  to  the  lines  of  constant  phase,  equal  to  zero: 


d[ut  —  <£(r)]  =  0  . 


(3.10) 


The  spatial  derivative  will  be  taken  with  respect  to  the  distance  £  in  the  direction 
h,  which  is  normal  to  the  wave  front.  (The  direction  n  varies  with  the  position 
of  the  wave  front.)  Equation  (3.10)  can  then  be  expanded  using  the  chain  rule: 


-  <j>{£))dt  +  -  <j>{i))dt  = 

dd> 

udt  +  —d£  =  0 


(3.11) 


The  propagation  speed  of  the  constant  phase  surface,  normal  to  the  surface,  is 
defined  as  cpd=^.  The  derivative  of  the  phase  with  respect  to  the  distance 
normal  to  the  constant  phase  ,  §|,  is  the  amplitude  of  the  phase  gradient,  |V<£|. 
Thus  the  phase  speed  cp  can  be  calculated  from  Equation  (3.11)  to  be 


Using  the  equation  for  the  active  intensity  and  potential  energy  \V<f>\  can  be 


expressed  as 


(3.13) 


where  V  —  jjr-P2  is  the  potential  energy.  Thus,  a  more  convenient  form  for  c. 


cp  —  2c  _  , 


which  expresses  cp  in  terms  of  two  measurable  quantities. 


(3.14) 


It  should  be  noted  that  cp  is  dependent  on  position.  Consequently,  at  each 
point  in  space  and  along  any  wave  front  the  phase  speed  may  vary.  Further, 
the  phase  speed  is  called  resultant,  to  make  a  clear  distinction  between  two  view 
points  of  an  acoustic  field.  Huygen’s  principle  concludes  that  an  acoustic  field 
can  be  represented  as  a  result  of  the  interference  of  waves  which  travel  at  speed 
c.  Thus  the  resultant  wave  propagates  at  the  resultant  phase  speed  not  at  the 
speed  c  at  which  the  interfering  waves  travel. 

A  quantity  often  referred  to  is  a  wave  number.  For  plane  and  spherical 
waves  the  wave  number  is  j.  In  general  we  can  discuss  a  resultant  wave  number 
kp  =  t-,  which  varies  in  space.  The  concept  of  a  spatially  dependent  wave  number 
is  not  easily  accepted  because  of  a  common  mathematical  tool  of  transforming 
real  space  into  wave  number  space.  The  transform  spaces  have  a  one  to  one 
mapping.  But  the  wave  number,  one  of  the  transform  variables,  is  independent 
of  real  space,  the  other  transform  variable. 

From  the  viewpoint  of  a  wave  number  decomposition  of  an  acoustic  field, 
where  the  wave  number  and  space  is  a  transform  pair  such  that  a  quantity  with 
the  units  of  wave  number  is  an  independent  variable,  the  variation  of  the  pressure 
phase  is  a  result  of  the  superposition  of  wave  numbers.  The  amplitude  of  the 
contribution  from  each  wave  number,  k,  is  independent  of  space;  where  as  the 


vx.-xo.  ->  .>  wo,  _>  V.’  \  ■/ 


resultant  wave  number  kp  =  f-,  which  results  from  a  superposition  of  all  the 
wave  numbers  in  the  field  and  is  a  directly  measurable  quantity  which  varies  in 
space.  Thus,  kp  defines  a  local  wave  number  in  contrast  to  k  which  is  a  global 
wave  number  that  is  averaged  over  all  space. 

It  is  often  convenient  to  split  a  field  into  a  superposition  of  waves,  wave 
numbers,  point  sources,  or  some  other  simplified  form.  However,  the  Helmholtz- 
Huygen  integral  states  that  for  a  given  pressure  and  velocity  distribution  there  is 
not  a  unique  source  distribution  which  creates  the  field.  Thus,  there  is  an  infinite 
set  of  realizable  decompositions  of  an  acoustic  field.  The  viewpoint  of  this  thesis 
is  to  visualize  an  acoustic  field  in  its  resultant  form,  without  decomposing  it  into 
different  forms. 

3.6.  Differential  Equations  for  P(r)  and  <j>(r)  in  Terms  of  the  Wave 
Properties 

The  resultant  wave  number  can  be  used  to  simplify  the  differential  Equations 
(3.5)  and  (3.6),  which  contains  the  term  V</>-  V<j>.  The  resultant  wave  number  in 
terms  of  the  phase  gradient  is 

kp  =  -  =  ~  =  \V<H  ■  (3.15) 

CP  |V0| 

In  vector  form  this  is 

kp  =  V<f>  .  (3.16) 

Substituting  this  into  Equation  (3.5)  and  using  the  fact  that  kp  •  kp  =  k%  results 
in  the  following  differential  equation  for  the  pressure  amplitude. 
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Equations  (3.17)  and  (3.18)  are  interesting  because  the  pressure  amplitude  at 
every  point  in  space  is  a  simple  function  of  the  wave  number,  k ,  and  the  resultant 
wave  number,  kp ,  at  that  point.  Consequently,  any  field,  even  a  complicated 
nearfield,  can  be  viewed  microscopically  on  a  very  simple  basis. 

Using  Equation  (3.17)  and  (3.18),  let  us  investigate  the  example  of  a  plane 
wave  where  the  phase  distribution  can  be  specified  as  <f>( x)  =  kx  which  are  waves 
propagating  at  speed  c.  Consequently,  kp  =  k.  Equation  (3.17)  and  (3.18) 


become 


rP  =  0  , 


2k— P  —  0  . 
ox 


(3.19) 

(3.20) 


The  only  solution  to  both  equations  is  P  =  constant,  which  is  the  form  for  a 
plane  wave. 

The  significance  of  Equation  (3.17)  and  (3.18)  is  that  for  a  specific  phase 
distribution  the  pressure  amplitude  is  specified  by  these  equations;  likewise, 
for  a  specific  pressure  amplitude  distribution  the  equations  specify  the  phase 
distribution.  Consequently,  if  the  potential  energy  is  known  as  a  function  of 
space  the  active  and  reactive  intensity  can  be  calculated  using  the  two  differential 
equations. 

3.7.  Example:  Plane  Wave  Incident  on  a  Rigid  Boundary 

Let  us  tie  together  some  of  the  ideas  on  the  wave  form  by  the  example  of  a 
plane  wave  incident  on  a  rigid  boundary,  Figure  3.1.  The  incident  pressure  is 


p.(x,y,t)  =  Ait 


—  A  .-jivt  —  k^xainO+yaine)) 


(3.21) 


The  total  pressure  is  the  sum  of  the  incident  wave  and  a  wave  reflected  from  the 
boundary.  The  velocity  can  be  calculated  from  the  pressure,  so  the  total  pressure 
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Figure  3.1  A  sketch  of  the  configuration  for  a  plane  wave  incident 
on  a  rigid  boundary.  The  phase  speed,  ep,  and  the 
group  velocity,  c  sm(0),  are  clearly  marked. 
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and  particle  velocity  are 


p(x,  y,  t)  =  2 Acos{kycosd)ej{ut-kxain^  , 
u{x,y,t)  =  i\^j^cos{kycos6)e^ut-kxain^} 

+  j[— —sin{kyCOs0)ei^t-*--kxaine'>\  . 
pc 

Using  these  two  equations  the  intensity  and  energy  densities  are 

_  2  A2 

I(x,y)=i - sin0cos2(2kycos6)  , 


(3.22) 


(3.23) 


(3.24) 


A2 

Q(x,y)  —  j — cosffsin(2kycosff)  , 
pc 

(3.25) 

A2 

V[x,y)  =  — -cos2  (2kycosff)  , 
pc* 

(3.26) 

A2  1 

T(x,y)  =  — -[-  —  cos2ffcos(2kycosff)\  . 

(3.27) 

There  are  some  features  of  these  equations  which  should  be  noted  and  are 
sketched  in  Figure  3.2.  The  active  intensity  is  parallel  to  the  rigid  boundary; 
consequently,  the  wave  fronts  are  surfaces  perpendicular  to  the  boundary 
propagating  in  the  positive  x  direction.  The  reactive  intensity  only  has  a 
component  perpendicular  to  the  boundary,  so  the  pressure  fluctuates  only  as 
a  function  of  the  distance  above  the  boundary.  This  can  also  be  visualized  as  a 
vertical  standing  wave. 

The  phase  speed  is  calculated  from  Equation  (3.14): 

cp(  i  =  2c2-=t  =  — —  .  3.28 

PK  ’  j/|  stuff  v  ' 

As  stated,  the  direction  of  the  resultant  wave  propagation  is  parallel  to  the 
boundary,  which  is  marked  Figure  3.2.  Two  extremes  can  be  looked  at.  First 
the  case  of  ff  =  \  which  is  the  grazing  angle  of  the  incident  wave  so  there  is  no 
reflected  wave.  The  wave  speed  for  ff  =  ^  is  c  Unce  there  are  no  interference 
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affects.  However  for  the  case  of  6  =  0,  cp  =  oo  which  is  the  case  of  the  incident 
wave  being  normal  to  the  boundary  so  there  is  no  component  horizontal  to  the 
boundary. 


For  this  example  the  wave  vector  is 


kp  =  i—sind  , 


(3.29) 


which  is  independent  of  space.  This  brings  up  an  interesting  point.  Going  back 
to  Equation  (3.18)  if  1  is  independent  of  space  V2d>  =  0,  then  Equation  (3.18) 


becomes 


2 kp  •  VP  =  0 


(3.30) 


Consequently,  the  active  intensity  and  reactive  intensity  are  perpendicular 
everywhere  when  the  resultant  wave  speed  is  independent  of  space. 

Another  issue  which  needs  to  be  addressed  is  the  speed  of  energy  propagation. 
Conventional  views  are  that  the  active  intensity  represents  the  energy  flow; 
however,  this  example  clearly  shows  a  problem  with  this  interpretation.  With  the 
incident  angle  6  the  component  of  the  incident  wave  which  propagates  parallel  to 
the  rigid  boundary  is  c  sinO.  The  implication  of  the  conventional  view  is  that  the 
energy  would  propagate  with  the  wave  fronts.  In  this  example  the  wave  fronts 
propagate  faster  than  the  speed  of  sound.  However,  the  acoustic  energy  cannot 
travel  faster  than  the  energy  which  creates  the  acoustic  field  which  propagates  at 
the  speed  of  sound.  Therefore,  in  this  example  the  group  velocity,  which  is  the 
average  speed  of  the  energy  propagation,  is  the  speed  of  the  energy  propagation 
parallel  to  the  rigid  boundary,  which  is  less  than  the  speed  of  sound. 
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Chapter  4 


CONCEPTS  INCORPORATING  ENERGY  TRANSFER 
AND  WAVE  PROPERTIES 


4.1.  Introduction 

The  previous  two  chapters  addressed  the  implications  of  the  intensity  vectors 
in  regard  to  energy  transfer  and  to  wave  motion.  Both  of  these  interpretations 
must  be  used  to  take  full  advantage  of  the  intensity  analysis.  This  chapter 
will  discuss  an  explanation  of  the  reactive  intensity  contribution  to  energy 
propagation  which  is  based  on  the  ideas  in  both  Chapters  2  and  3.  The  group 
velocity  and  practical  aspects  of  the  reactive  intensity  will  also  be  discussed. 
Finally,  the  example  of  the  active  intensity  vortex  is  explored,  from  the  viewpoint 
of  the  wave  properties  and  energy  transfer. 

4.2.  A  Comparison  of  Wave  Front  and  Energy  Propagation 

Fundamental  explanations  of  wave  phenomena  often  use  the  special  case  of 
the  plane  wave  to  show  that  energy  propagates  with  the  wave  fronts.  However, 
combining  the  interpretations  of  Chapters  3  and  4,  the  energy  travels  with  speed 
c  in  the  direction  of  the  instantaneous  intensity  which  in  a  general  wave  is  not 
in  the  direction  of  the  wave  fronts.  It  is  difficult  to  conceptualize  the  idea  that 
energy  does  not  propagate  with  the  waves,  because  intuitively,  in  a  pressure  wave, 
at  a  time  or  in  a  position  with  zero  pressure,  no  acoustic  energy  is  expected. 
However,  the  contribution  of  the  reactive  intensity  is  not  considered  in  such  an 
argument.  With  the  wave  properties  linked  to  the  active  intensity  in  Chapter  3 
a  better  physical  understanding  of  the  reactive  intensity  can  now  be  introduced 
to  clarify  this. 
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The  differential  equations  for  P[r )  and  4>{r) ,  Equations  3.5  and  3.6,  specify 
that  if  the  resultant  phase  speed  is  different  than  the  plane  wave  phase  speed, 
then  there  is  a  pressure  gradient.  The  pressure  gradient  means  that  there  is  a 
non-zero  kinetic  energy  component  due  to  Ur. 

Keeping  these  ideas  in  mind  let  us  consider  a  plane  wave.  The  kinetic  energy 
is  in  phase  with  pressure,  because  the  pressure  amplitude  does  not  vary  in  space, 
UT  =  0.  Figure  4.1  represents  the  time  dependent  potential  energy  and  kinetic 
energy  at  times  t  and  some  small  increment  of  time  At  later.  The  energy  is 
contained  in  the  regions  marked  (I),  (II),  (III),  and  (IV).  At  time  t  +  At  the 
energy  has  progressed  to  the  right  at  speed  c  in  the  packets  (I),  (II),  (III),  and 
(IV).  This  is  one  simple  way  to  demonstrate  that  the  energy  propagates  with  the 
plane  waves. 

Now  consider  the  case  where  the  resultant  waves  do  not  travel  at  speed  c 
(cp  ^  c),  so  that  there  is  kinetic  energy  due  to  the  components  Ua  and  UT.  Figure 
4.2  represents  the  time  dependent  energy  densities  as  a  function  of  distance  for 
times  t  and  t  +  At.  The  reference  points  which  we  call  wave  fronts,  propagate 
the  distance  cpAt,  thus  the  packets  between  pressure  zeros  labeled  (I),  (II),  and 
(III)  propagate  a  distance  cvAt.  The  energy  propagates  at  speed  c  which  in  this 
case  is  not  with  the  packets  thus  there  is  more  energy  at  point  (1)  in  the  kinetic 
energy  associated  with  Ur  than  in  the  wave  form.  In  order  to  maintain  the  shape 
of  the  propagating  waves  (I),  (II),  and  (III),  kinetic  energy  associated  with  UT  at 
point  (2)  is  shifted  into  Ua  and  p. 

To  summarize  this  argument,  it  must  be  remembered  that  for  the  plane  wave 
the  energy  travels  with  the  waves  because  there  is  no  energy  associated  with  Ur. 
However,  with  a  non-zero  component  of  the  kinetic  energy  due  to  Ur  and  with 
cp  ^  c,  the  energy  propagates  by  shifting  between  the  kinetic  energy  due  to  Ua 
and  UT  so  the  wave  fronts  can  propagate  at  a  speed  other  than  the  speed  of 
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Figure  4.1  The  magnitude  of  the  pressure  and  velocity  as  a 
function  of  distance  in  the  direction  of  a  propagating 
plane  wave  at  times  t  and  t  +  At.  (At  =  ~)  The 
wave  packets  I,  II,  III,  and  IV  are  labeled  so  that  their 
progress  can  be  followed. 
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Figure  4.2  The  magnitude  of  the  pressure  and  the  velocities  Ua 
and  Ur  as  a  function  of  distance  in  the  direction  of 
a  propagating  wave,  where  the  pressure  amplitude 
decreases  in  the  propagation  direction  and  cp  7^  c, 
at  times  t  and  t  +  At.  (At  =  j)  The  wave  packets 
I,  II,  and  III  are  labeled  so  that  their  progress  can  be 
followed. 


energy  propagation,  c.  The  wave  fronts  are  points  of  reference  on  a  wave  form. 
Only  in  special  cases  where  Ur  =  0  is  the  propagation  characteristics  of  the  wave 
front  the  same  as  the  energy  propagation. 

4.2.1.  Group  Velocity 

The  group  velocity  was  devised  to  quantify  the  speed  of  energy  propagation 
through  a  region  of  space.  This  is  a  deceiving  concept  because  it  has  already  been 
stated  that  energy  propagates  at  speed  c  in  the  direction  of  the  instantaneous 
intensity.  However,  we  have  seen  with  the  instantaneous  energy  paths  that  energy 
does  not  propagate  over  a  single  path  for  all  times  during  one  period.  Thus  a 
global  description  of  energy  propagation  is  in  order. 

The  acoustic  wave  guide  is  the  best  example  of  the  intent  of  group  velocity. 
In  an  acoustic  wave  guide;  i.e.,  a  duct,  the  propagation  speed  of  wave  fronts, 
cp,  can  be  greater  than  the  speed  of  sound.  However,  the  active  intensity  does 
not  include  all  of  the  energy  flux  in  the  duct  since  there  is  reactive  intensity. 
Incorporating  the  energy  propagation  described  by  the  reactive  intensity  results 
in  the  speed  of  energy  propagation  down  the  duct  which  is  less  than  c.  Thus 
the  group  velocity  is  defined  in  terms  of  the  speed  at  which  energy  at  one  cross 
section  of  the  duct  propagates  to  another  cross  section  down  the  duct. 

There  are  several  derivations  of  the  group  velocity  (Lamb,  1904;  Biot, 
1957;  Whitham,  1961;  Jackson,  1975).  However,  a  basic  assumption  of  all  the 
derivations  specifies  the  farfield  and  time  averages  for  mathematical  ease.  The 
conceptual  issue  is  the  validity  of  describing  a  time  varying  process  of  energy 
propagation  with  a  time  independent  quantity.  Presently,  we  are  only  able 
to  quantify  the  effect  of  the  time  dependent  process  of  energy  propagation, 
but  cannot  conveniently  or  meaningfully  describe  the  time  dependent  energy 
propagation  itself. 


The  actual  path  of  energy  propagation  are  themselves  open  for  interpretation 
and  in  fact  there  is  not  a  single  path,  Section  2.5.  Consequently,  investigations  of 
energy  flow  with  the  instantaneous  intensity  is  not  practical  and  may  lack  precise 
meaning.  The  group  velocity  is  a  compromise  and  is  a  means  to  quantify  the 
idea  of  energy  propagation;  however,  direction  cannot  be  attached  to  the  group 
velocity  as  we  cannot  attach  direction  of  the  energy  propagation  with  the  active 
intensity  vector. 

4.3.  Reactive  Intensity:  A  Practical  Viewpoint 

The  reactive  intensity  has  received  little  attention  in  the  literature  and  is 
generally  ignored  since  its  contribution  to  the  time  average  of  the  instantaneous 
intensity  is  zero  at  a  given  point.  However,  as  emphasized  in  Chapters  3  and  4 
the  vector  does  have  physical  significance  and  thus  helps  characterize  the  acoustic 
field. 

One  refreshing  article  was  published  by  Wolfgang  Schmidt  (1985),  in  which  he 
used  the  intensity  analysis  in  the  nearfield  of  rotating  machinery.  The  interesting 
aspect  of  his  analysis  was  that  he  recognized  that  the  component  of  the  particle 
velocity  in  phase  quadrature  with  the  pressure  is  dominant  in  the  nearfield  of  an 
acoustic  source.  Thus  by  electronically  shifting  the  measured  velocity  signal  90 
degrees  Schmidt  was  able  to  measure  the  reactive  intensity  to  characterize  faulty 
equipment. 

Indeed  in  the  nearfield  the  reactive  intensity  dominates  the  active  intensity, 
often  by  10  to  15  dB.  Thus  a  significant  part  of  the  acoustic  energy  in  the 
nearfield  is  described  by  the  reactive  intensity.  Let  us  revert  for  a  minute  to  a 
pure  standing  wave  with  a  maximum  SPL=  140  dB.  The  active  intensity  will 
be  zero  so  there  is  no  net  energy  flux,  which  in  a  viewpoint  of  noise  control  is 
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a  desirable  result;  however,  a  person  in  such  an  acoustic  field  will  risk  serious 
hearing  loss  if  exposed  for  substantial  periods  of  time. 

Consequently,  if  we  wish  to  take  a  noise  control  engineering  approach  to  a 
noise  source  we  can  say  that  we  are  interested  in  the  sound  which  a  person  is 
exposed  to,  which  we  need  to  somehow  quantify.  If  the  person  is  in  the  farfield  of 
a  source,  then  the  power  output  of  a  source  is  a  valid  description  of  the  acoustic 
energy  to  which  the  individual  will  be  exposed.  However,  if  an  individual  is 
situated  in  the  nearfield  of  the  acoustic  source  then  the  energy  which  is  described 
by  the  reactive  intensity  can  be  significant  if  not  dominant.  Thus  quantifying  the 
source  radiation  by  the  active  intensity  alone  is  an  inadequate  characterization 
of  the  energy  to  which  an  individual  is  exposed  in  the  nearfield  of  a  source. 

4.4.  Active  Intensity  Vortex 

The  active  intensity  vortex  has  had  a  dubious  rule  in  the  development  of 
intensity  over  the  past  20  years.  The  active  intensity  vortex  is  a  region  where  the 
active  intensity  vectors  form  a  vortex  pattern  which  is  commonly  seen  in  fluid 
flow.  In  1964,  Schultz,  Smith  and  Malme  made  some  measurements  in  which 
they  concluded  that  the  active  intensity  was  forming  circular  paths.  In  a  paper 
published  in  June  1975  where  they  revealed  this  conclusion,  they  frankly  said 
that  these  were  results  they  “did  not  have  the  nerve  to  present  at  that  time, 
because  they  did  not  make  sense.”  Since  then  there  have  been  numerous  authors 
who  have  presented  evidence  and  explanations  for  the  vortex  (Kristiansen  1981; 
Elko  1984;  Tichy  and  Mann  1985;  Waterhouse,  Crighton,  and  Ffowcs-Williams 
1987;  Waterhouse  1987;  Mann,  Tichy,  and  Romano,  1987). 

Understanding  the  vortex  has  been  a  pivotal  issue  in  much  of  the  work 
presented  in  this  thesis.  The  vortex  developed  into  an  interesting  conversation 
piece  because  it  presented  a  strong  case  against  conventional  views  of  the  active 


intensity  representing  total  power  flow.  Such  a  perspective  says  that  the  energy 
in  the  vortex  is  circulating  around  the  center  point,  thus  representing  energy 
which  is  locally  trapped  in  a  region  of  space.  This  concept  was  indeed  puzzling 
and  required  an  understanding  of  the  intensity  vector  which  was  not  readily 
available  in  the  literature, 

4.4,1.  The  Conditions  for  the  Vortex  Formation 

The  existence  of  the  vortex  pattern  comes  about  because  the  curl  of  the 
active  intensity  vector  is  non-zero.  This  allows  the  pattern  to  be  formed  by  the  I 
vector.  Before  addressing  the  issue  of  energy  flux,  let  us  use  the  wave  structure 
interpretation  detailed  in  Chapter  3  to  explain  the  vortex. 

The  configuration  of  three  point  sources  will  be  used  for  computational  ease. 
The  sources  are  in  a  straight  line  with  18  cm  spacing  and  with  relative  amplitudes 
of  +1,-1  and  +1.  The  randomly  chosen,  frequency  of  2070  Hz  is  used.  Figure  4.3 
is  the  active  intensity  superimposed  with  wave  fronts  of  equal  phase  separation. 
There  are  two  regions  where  the  wave  fronts  intersect.  The  left  region  is  enlarged 
in  Figure  4.4  along  with  the  reactive  intensity.  This  region  will  be  explored  in 
detail. 

There  are  two  points  in  the  vortex  region  that  are  of  special  significance:  the 
vortex  center  and  the  saddle  point.  The  reactive  intensity  points  into  the  center 
of  the  vortex,  meaning  there  is  a  pressure  minimum  at  the  center.  In  Figure  4.5 
various  quantities  are  plotted  along  a  line  passing  through  the  saddle  point  and 
vortex  center.  The  potential  energy,  active  intensity,  and  reactive  intensity  are 
zero  at  the  vortex  center.  However  neither  the  curl  of  I  or  the  kinetic  energy  are 
zero  or  a  maximum  at  the  center. 
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Figure  4.3  The  calculated  active  intensity  and  wave  fronts  above 
three  point  sources  separated  by  18cm  at  2070Hz  with 
the  center  source  out  of  phase.  The  vortex  region  is 
circled. 


Il 

p, 


^  \  *\  ^  « —  * —  *•""  t-*'  * —  « —  ** —  *\  \ 

\  \  \  ^«—  r'  S  S'  S  ^  +~  *^\\ 

1  1  \  \-  //  //  /-\U 

r  r  i  Tfi  i  i  i  i  j  '  \  i  i 
] k//W'-\\\\\4/;r 
/  / 
/  /  ^  ✓*  / 


K#* 
•JtJjV, 


I®® 


i  tude 


i.ri.i  i, 


- cur  I  of  intens i ty 

— .  potential  energy 

-  active  intensity 

.  \\ 

\ 


-  kinetic  energy 

-  reactive  intensity 


\ - 

\  \ 


\/ 

/  \ 

\  /  \ 


saddle  point ' 


d  i  s  t  once  cm  vortex  center 


Figure  4.5  The  scaled  magnitude  of  active  intensity,  reactive 
intensity,  potential  energy,  kinetic  energy,  and  curl  of 
the  active  intensity  along  a  line  through  the  saddle 
point  and  vortex  center. 
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At  the  saddle  point  the  curl  of  I  is  zero.  The  saddle  point  is  always  present 
in  a  vortex  region  because  it  is  the  point  where  the  circular  pattern  of  the  vortex 
is  isolated  from  the  rest  of  the  acoustic  field. 

The  most  interesting  quantities  are  the  particle  velocity  and  pressure  phase. 
As  already  stated,  the  total  kinetic  energy  has  no  significant  properties  in  the 
vortex  region;  however,  the  two  components  Ua  and  Ur  provide  an  interesting 
clue.  At  the  saddle  point  Ua  =  0  since  7  =  0.  Figure  4.6  is  a  close  up  of  JJa  and  Ur 
around  the  vortex  center.  The  relative  magnitude  normalized  to  the  maximum 
of  each  quantity  is  plotted  so  that  both  quantities  could  be  represented  together. 
At  the  vortex  center  Ur  has  a  relative  minimum  and  Ua  has  a  relative  maximum. 
The  peak  in  Ua  has  special  significance.  Since  there  cannot  be  an  infinite  amount 
of  energy  at  a  point,  Ua  has  a  finite  limit  at  the  center,  call  it  /c.  Now  solve  for 
V4>  in  the  equation  for  Ua,  Equation  (2.6): 

V(£  =  =  up—  .  (4.1) 

However  at  the  vortex  center  P  —  0;  consequently  at  the  vortex  center 
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The  spatial  phase  gradient  is  plotted  in  Figure  4.7.  The  resolution  of  the  plot  is 
sufficient  to  see  that  Vcf>  is  going  to  infinity.  This  means  that  the  pressure  phase  is 
discontinuous  at  the  vortex  center.  However,  it  is  true  that  p(r,f)  is  smooth  and 
continuous  through  the  vortex  center,  goes  to  zero  at  the  discontinuity;  and  is 
positive  on  one  side  of  the  zero  and  negative  on  the  other.  This  is  also  observed  in 
plane  waves  at  nodal  points  where  the  physical  quantities  p  and  u  are  continuous 
but  the  phase  is  discontinuous  by  |. 

The  phase  discontinuity  explains  the  vortex  pattern  and  represents  the 
condition  for  the  formation  of  a  vortex.  In  general,  a  phase  discontinuity  occurs 
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at  the  intersection  of  different  surfaces  of  constant  phase.  Such  intersections  can 
occur  along  a  line,  or  at  a  point,  so  that  the  phase  discontinuities,  and  thereby 
the  vorticies,  must  be  on  a  line  or  at  a  point. 

There  are  many  examples  of  phase  discontinuities.  Let  us  first  consider  a 
dipole.  On  the  plane  between  the  two  sources  the  pressure  is  zero  and  the  phase 
is  discontinuous.  The  surfaces  of  constant  phase  must  be  perpendicular  to  this 
plane  of  zero  pressure  which  means  that  the  active  intensity  must  be  parallel  to 
the  dividing  plane.  Consequently,  a  phase  discontinuity  which  forms  a  plane  will 
not  produce  a  vortex  pattern. 

Now  consider  the  case  when  the  pressure  phase  is  discontinuous  on  a  line. 
(This  is  also  a  pressure  null  line.)  The  different  surfaces  of  constant  phase  will 
extend  radially  from  the  pressure  null  line.  The  same  will  occur  if  the  phase 
is  discontinuous  at  a  point.  Since  the  active  intensity  is  perpendicular  to  the 
wave  fronts  the  active  intensity  will  form  a  circular  pattern,  a  vortex,  around 
the  line  or  point  of  phase  discontinuity.  Clearly,  the  vortex  physically  represents 
a  resultant  wave  travelling  around  a  phase  discontinuity  only  when  the  phase 
discontinuity  (also  where  the  pressure  is  zero)  is  a  line  or  a  point. 

This  description  of  the  vortex  is  consistent  with  the  wave  front  interpretation 
in  Chapter  3  and  is  in  essence  the  observation  made  by  Kristiansen  (1981).  The 
discontinuity  could  have  been  implied  by  constructing  wave  fronts  from  the  active 
intensity.  Connected  with  the  wave  fronts  is  the  phase  speed  which  is  presented  in 
Figure  4.8.  This  shows  that  the  speed  of  the  wave  fronts  increases  with  distance 
from  the  vortex  center.  Also  there  are  regions  where  the  waves  progress  slowly 
and  others  where  they  propagate  very  quickly. 


4.4.2.  Energy  Flux  in  the  Vortex 

Having  resolved  the  meaning  of  the  vortex  formation,  the  more  significant 
issue  is  the  energy  flux  in  the  vortex.  Statements  have  been  made  that  the 
energy  in  the  vortex  arrives  during  the  transient  build  up  of  the  vortex,  and 
thus  flows  in  closed  paths  until  the  sound  sources  are  turned  off  and  the  field 
decays  (Kristiansen  1981;  Waterhouse,  Crighton,  and  Ffowcs- Williams  1987; 
Waterhouse  1987).  Such  interpretations  are  in  direct  conflict  with  the  concepts 
of  energy  transfer  presented  in  Chapter  2. 

The  transient  build  up  of  the  vortex  can  be  investigated  by  considering 
the  time  dependent  pressure  and  particle  velocity.  The  transient  response  is 
incorporated  by  considering  the  contribution  of  each  point  source  only  after  the 
first  wave  front  from  each  source  has  reached  the  point  where  the  pressure  and 
velocity  are  being  calculated  (Junger,  1966).  This  is  expressed  by  including  the 
unit  step  function  in  the  summation  such  that  from  the  time  the  sources  are 
turned  on  until  the  wave  traveling  at  speed  c  reaches  a  point  the  contribution  to 
the  pressure  and  velocity  is  zero: 

Pir,t)  =  Px{r,t)H[ct  -  ri)  +  p2{r,t)H(ct  -  r2)  +  -  r3)  ,  (4.3) 

u(r,t)  =  ux{r,t)H{ct  —  ri)  +  u2(r,t)H(ct  —  r2)  +  u3(r,t)H(ct  -  r3)  .  (4.4) 

The  instantaneous  intensity  can  subsequently  be  calculated.  Let  us  look  at  the 
part  of  the  instantaneous  intensity  containing  the  active  intensity,  7(r)  cos2{ut  — 
4>).  Figure  4.9  contains  plots  of  this  term  for  equal  increments  of  time. 

In  Figure  4.9(a)  the  waves  from  the  two  closest  sources  have  already  arrived 
in  the  region  of  the  vortex.  At  the  lower  right  corner  of  Figure  4.9(b)  the 
contribution  from  the  third  source  is  seen.  The  following  plots  show  the 
progression  of  the  wave  through  the  region.  It  is  important  to  notice  that  once 
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Figure  4.9  The  contribution  of  the  active  intensity  to  the  instan¬ 
taneous  intensity  during  the  vortex  buildup.  The  time 
increment  is  0.0387. 


the  third  wave  arrives  at  a  point  the  orientation  of  the  active  intensity  vector  is 
in  the  pattern  of  the  vortex.  There  is  no  transient  trapping  of  the  energy  in  the 
vortex. 

A  more  convincing  argument  is  to  view  the  instantaneous  intensity  in  the 
steady  state  condition.  Figure  4.10  is  a  sequence  of  snapshots  over  one  cycle  of 
the  instantaneous  intensity  in  the  vortex  region  at  equal  time  increments.  These 
plots  bear  no  resemblance  to  the  vortex  pattern,  because  the  reactive  intensity 
has  a  significant  contribution  to  the  energy  transport. 

The  transfer  of  energy  can  also  be  described  by  the  instantaneous  flux  paths 
which  were  described  in  Chapter  2.  These  are  drawn  in  Figure  4.11.  Figure 
4.11(a)  is  the  region  of  the  vortex  and  Figure  4.11(b)  is  an  expanded  view  which 
follows  paths  for  more  than  one  cycle.  The  seven  starting  positions  are  circled, 
and  two  times  in  one  period  are  chosen  to  follow  the  energy.  Figure  4.11(a) 
clearly  shows  paths  which  do  not  resemble  the  vortex,  in  fact  all  the  paths  leave 
the  plotted  region.  Figure  4.11(b)  shows  that  the  paths  leaving  the  vortex  region 
continue  and  never  return.  Consequently  energy  is  not  trapped  in  the  vortex. 

Another  argument  against  the  belief  of  energy  being  trapped  in  the  vortex 
was  made  by  Romano  (1986).  Romano’s  argument  is  based  on  the  interpretation 
of  the  acoustic  field  being  the  superposition  of  energy  contribution  from  in  this 
case  the  three  sources.  Following  through  this  argument  in  turn  concludes  that 
energy  is  not  trapped  in  the  vortex. 
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Figure  4.10  Snapshots  of  the  instantaneous  intensity  in  the  cortex 
region  for  one  complete  cycle  of  the  energy  (y).  The 
time  increment  is 


Chapter  5 


EXPERIMENTAL  PROCEDURE  AND 
DATA  PROCESSING 


5.1.  Introduction 

This  chapter  discusses  the  experimental  facilities  and  the  processing  of  the 
measured  data.  The  emphasis  of  this  thesis  is  understanding  and  interpreting 
acoustic  intensity;  consequently,  there  was  little  development  of  the  measurement 
procedures  beyond  the  available  instrumentation. 

5.2.  Experimental  Facilities 

The  experimental  facilities  used  were  developed  by  Gary  Elko  and  described 
in  his  dissertation  (1984).  The  system  is  controlled  by  a  PDPll/34  computer, 
fully  automating  the  measurement  procedure.  The  data  are  transferred  to  a  VAX 
780  computer  for  phase  and  impedance  calculations,  plotting,  and  permanent 
storage  on  magnetic  tape. 

The  measurements  are  performed  in  a  2.36m  by  1.51m  by  1.89m  room  with 
hard  walls.  After  installing  foam  wedges  to  make  the  room  anechoic  the  available 
uninterrupted  interior  space  is  1.4m  by  0.85m  by  1.6m.  The  wedges  provide 
an  anechoic  environment  to  the  lowest  frequency  of  300Hz,  which  establishes 
the  lower  limit  for  measurements  requiring  an  anechoic  environment.  A  three- 
dimensional  motion  system  installed  in  the  room  moves  the  intensity  probe  to 
any  position  within  the  room.  The  motion  system  is  controlled  by  the  PDPll/34 
computer  through  a  RS232  serial  communication  line. 

The  intensity  probe  consisted  of  four  B&K  ^  inch  pressure  microphones.  The 
four  microphones  are  grouped  in  two  pairs.  Each  pair  of  microphones  measure 
the  velocity  in  a  given  direction.  The  two  pairs  are  orientated  so  that  they 


measure  perpendicular  components  of  the  particle  velocity.  The  minimum  center 
to  center  spacing  of  the  microphones  is  17mm.  This  establishes  the  upper  limit 
of  the  measurement  system  at  4500Hz  for  the  microphone  spacing  to  be  within 
about  2  percent  of  the  acoustic  wavelength. 

The  two  microphone  pairs  were  phase  calibrated  using  the  modified  switching 
technique  in  a  free  field  as  outlined  by  Elko  (1984).  This  gave  a  corrected  phase 
with  the  accuracy  less  than  0.1  degrees,  which  is  a  conservative  estimate.  The 
phase  shift  introduced  by  the  microphone  cartridge  was  compensated  for  in  the 
calculation  of  the  cross-spectrum.  The  shift  was  a  constant  time  delay  so  the 
frequency  shift  was  linear.  The  calibration  curves  were  fitted  with  a  least  squares 
algorithm  to  the  form  Af  -f  B,  making  it  easy  in  the  data  processing  to  calculate 
the  phase  mismatch  at  each  frequency.  An  example  of  the  phase  calibration  is 
given  in  Figure  5.1. 

The  microphone  signals  are  amplified  50dB  to  70dB  so  that  the  signals 
are  well  within  the  dynamic  range  of  the  A/D  converter.  The  four  pressure 
signals  are  then  passed  through  anti-aliasing  filters  with  a  48dB  per  octave  roll 
off.  The  amplified  and  filtered  signals  are  input  to  the  A/D  converter  in  the 
PDPll/34  computer.  The  PDP  digitally  sample  the  four  signals  and  using  an 
array  processor  calculates  the  1024  point  FFT’s  and  the  512  point  auto-spectra 
and  cross-spectra.  The  intensity  vectors,  energy  densities,  and  the  bias  errors  are 
calculated  from  the  spectra  of  the  two  microphone  pairs.  The  spectrum  analysis 
is  detailed  in  Elko’s  dissertation. 

The  PDP  fully  automates  the  measurement  process  by  moving  the  micro¬ 
phone  probe  to  specified  positions  in  the  chamber  and  acquiring  data.  The 
process  takes  about  three  and  a  half  hours  for  300  data  points  when  a  single 
frequency  is  being  used.  Often  measurements  were  made  around  an  obstacle  in 
the  scanning  plane.  For  such  cases  the  measurement  was  performed  in  separate 
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segments  which  followed  the  contour  of  the  obstacle.  The  data  files  from  the 
segments  were  subsequently  combined  so  that  the  the  entire  field  around  the 
obstacle  was  represented. 

Measurements  were  performed  at  single  frequencies.  To  speed  up  some 
measurements  multiple  frequencies  were  input  to  the  source;  however,  at 
frequencies  far  apart.  Broad  band  signals  were  avoided  because  of  the  long 
averaging  time  needed  for  adequate  accuracy.  Broad  band  signals  also  produce 
difficulties  because  there  are  low  signal  levels  at  desired  frequencies  bands  since 
the  energy  is  spread  over  a  large  frequency  range. 

Many  measurements  can  be  performed  on  a  symmetric  plane  so  that  the 
intensity  components  are  only  within  the  plane.  For  these  cases  the  two 
dimensional  four  microphone  probe  is  ideal.  When  the  intensity  is  three 
dimensional  the  measurements  are  made  with  the  probe  in  one  orientation,  first 
measuring,  for  instance,  the  x  and  y  components,  then  the  probe  is  rotated  90 
degrees,  now  measuring,  for  instance,  the  y  and  z  components.  The  data  are 
then  combined  and  the  common  measured  component  is  compared  to  monitor 
the  accuracy  of  the  measurement.  The  two  probe  orientations  are  measured  in 
separate  scans,  so  it  is  crucial  to  guarantee  that  the  measurement  positions  are 
accurately  repeated.  There  have  been  some  probe  configurations  developed  to 
measure  three  components  simultaneously;  however,  these  are  prone  to  diffraction 
errors  and  for  mechanical  reasons  the  scan  cannot  be  made  close  to  surfaces.  The 
method  of  two  probe  orientations  has  been  accurate  and  measurements  can  be 
made  within  2.5cm  of  a  surface. 

5.3.  Plotting  Symbols 

The  representation  of  the  intensity  and  energy  quantities  is  consistent  with 
conventional  representations.  The  active  and  reactive  intensity  are  represented 


by  arrows  pointing  in  the  direction  of  the  intensity  vector.  The  potential  energy  is 
represented  by  circles.  The  kinetic  energy  is  represented  by  parallelograms  which 
display  the  two  spatial  components  of  the  kinetic  energy  by  varying  the  height 
and  width  of  the  parallelogram.  The  center  of  each  symbol  is  centered  at  the 
measurement  points.  The  intensity  vectors  are  often  drawn  with  the  base  at  the 
measurement  point,  but  we  have  found  centering  the  vector  at  the  measurement 
point  provides  the  best  graphic  dynamic  range.  The  symbols  are  scaled  in  dB  to 
the  maximum  value  in  the  plotted  area  with  a  30dB  dynamic  plotting  range.  The 
maximum  value  is  indicated  on  the  top  of  the  plot.  The  dB  value  of  each  quantity 
is  given  as  a  numeric  value  representing  the  dB  level  below  the  maximum  value. 

Measurements  are  also  presented  for  three-dimensional  regions.  In  these 
plots,  the  in  plane  components  of  the  intensity  vector  are  represented  by  the 
direction  of  the  arrow,  and  the  out-of-plane  component  is  represented  by  the 
thickness  of  each  arrow.  The  thickness  is  scaled  to  |  the  actual  out  of  plane 
value.  The  arrow  is  solid  when  pointing  out  of  the  paper  and  empty  when 
pointing  into  the  paper. 

5.4.  Data  Processing 

During  the  measurements  the  active  intensity,  reactive  intensity,  potential 
energy,  kinetic  energy  and  bias  errors  are  calculated.  The  phase  distribution, 
resultant  phase  speed,  and  specific  acoustic  impedance  are  calculated  once  the 
data  are  transferred  to  the  VAX.  The  calculation  of  the  resultant  phase  speed 
and  specific  acoustic  impedance  is  straight  forward;  however,  the  pressure  phase 
calculation  will  be  explained  here. 

5.4.1.  Calculation  of  the  Pressure  Phase 

As  shown  in  Chapter  2  the  active  intensity  is  proportional  to  the  pressure 
squared  times  the  gradient  of  the  pressure  phase.  Thus  the  gradient  of  the 


pressure  phase  can  be  calculated  as  follows  from  the  active  intensity  and  potential 


energy: 
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From  the  gradient  over  a  grid  of  points  numerical  integration  can  be  used  to 
calculate  the  pressure  phase  (Imaich  and  Ohmi  1983).  Surfaces  of  constant  phase 
can  then  be  plotted  at  equal  phase  increments  so  the  wave  front  shapes  and 
propagation  characteristics  can  be  effectively  visualized. 

The  gradient  is  spatially  integrated  using  the  trapezoidal  rule.  A  schematic 
of  the  procedure  is  presented  in  Figure  5.2.  The  integration  begins  at  the  center 
to  minimize  accumulative  errors.  Lines  extending  vertically  and  horizontally 
from  the  center,  labeled  I,  II,  III,  and  IV,  are  the  first  integration  paths.  From 
these  four  lines  the  other  points  are  approached  along  two  paths,  as  indicated  by 
two  arrows  pointing  to  each  point.  The  two  valves  obtained  from  each  path  are 
averaged  for  the  phase  at  that  point.  The  procedure  continues  for  all  the  points 
resulting  in  a  phase  distribution  referenced  to  zero  at  the  center  of  the  grid. 

Measurement  areas  containing  obstacles,  require  segmentation  of  the  inte¬ 
gration  region.  The  areas  must  be  segmented  such  that  no  point  is  reached  by 
integrating  through  an  unmeasured  area.  This  also  means  shifting  the  origin  of 
integration.  The  sections  are  then  integrated  separately,  but  referenced  to  one 
section  so  the  phase  is  continuous  through  the  entire  plot.  Figure  5.3  gives  an 
exampie  of  the  segmentation,  origin  relocation,  and  connection  scheme. 

The  phase  resulting  from  this  process  is  in  radians  referenced  to  some 
arbitrary  zero  phase,  which  is  unrelated  to  the  source  phase.  However,  if  desired 
the  data  can  be  referenced  to  the  source  by  a  quick  measurement  of  the  phase 
between  the  pressure  at  the  source  and  one  point  in  the  measurement  region. 
Since  the  data  is  in  radians,  the  wave  fronts,  lines  of  constant  phase,  can  be 


Figure  5.2 


A  schematic  of  the  integration  paths  used  to  recon¬ 
struct  the  pressure  phase  from  the  measured  data. 


Figure  5.3 


A  schematic  of  the  segmentation  (labeled  1,2,3),  origin 
relocation  (circled),  and  connection  paths  (arrows)  for 
the  reconstruction  of  the  pressure  phase  around  an 
obstacle  in  a  measurement  plane. 


drawn  a  fraction  of  cycle  apart.  This  is  convenient  since  the  resultant  local 
wavelength  can  be  visually  related  to  the  free  field  wavelength  (A  =  y) 


5.4.2.  Specific  Acoustic  Impedance 


The  specific  acoustic  impedance  is  intrinsic  in  concepts  of  acoustic  fields.  It 
is  a  complex  function  defined  by  the  ratio  of  the  pressure  and  particle  velocity. 
Lahti  (1985)  develops  the  impedance  directly  in  terms  of  the  auto-spectra  and 
cross-spectrum  of  two  microphone  signals.  However,  it  is  convenient  to  express 
the  impedance  in  terms  of  the  active  and  reactive  intensity  since  the  spectra  may 
not  be  directly  available.  The  impedance  is  a  scalar,  but  it  can  be  calculated  from 
the  particle  velocity  in  a  given  direction.  For  instance,  it  is  usually  desired  to 
calculate  the  impedance  normal  to  a  particular  radiating  surface.  Consequently, 
the  complex  specific  acoustic  impedance,  Z_n,  in  the  direction  of  the  unit  vector 
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where  Tn  is  the  kinetic  energy  due  to  the  velocity  in  the  direction  h.  The  specific 
acoustic  impedance  can  be  measured  at  any  point;  however,  the  impedance  has 
significant  meaning  at  boundaries,  because  of  sources  and  sinks  of  energy.  Thus 
measuring  the  impedance  around  acoustic  sources  and  diffraction  boundaries  is 
a  useful  analysis  tool.  The  measurements  presented  in  Chapters  6  and  7  will 
clearly  show  some  examples  of  the  specific  impedance  measured  at  boundaries. 

There  are  some  examples  in  the  literature  of  using  the  specific  acoustic 
impedance.  They  are  geared  toward  general  absorption  characteristics  of 
materials.  However  in  this  thesis  a  point  by  point  analysis  of  the  specific  acoustic 
impedance  will  be  made. 


EXPERIMENTAL  ANALYSIS  OF  SOURCE-RESONATOR 


INTERACTIONS  AND  A  SIMPLE  SOURCE 


6.1.  Introduction 


This  chapter  focuses  on  analysis  of  source  and  resonator  interactions  using 
the  acoustic  intensity  measurement  technique.  Machines,  other  noise  producing 
mechanisms,  and  transducers  which  are  designed  to  radiate  sound  can  be 
classified  as  a  combination  of  sources  with  significant  structural  vibration  and 
cavities  which  resonate  and  absorb  energy.  For  much  of  the  analysis,  the  active 
and  reactive  intensity  is  qualitatively  used  to  visualize  the  wave  front  propagation 
and  pressure  amplitude  distribution.  The  real  and  imaginary  parts  of  the  specific 
impedance  show  in  detail  the  sources  and  absorbers  and  distinguishes  between 
them. 

This  chapter  also  compares  the  measured  radiation  impedance  and  radiated 
power  of  a  simple  source  as  a  function  of  the  source  height  above  a  rigid  boundary. 
These  measurements  show  that  measurement  errors  are  introduced  by  signal 
processing  equipment  inaccuracies  besides  phase  the  mismatch  of  the  signal 
processing  channels. 

6.2.  Two  In  Phase  Piston  Sources  Surrounding  a  Resonating  Tube 

The  first  example  is  a  configuration  sketched  in  Figure  6.1.  There  are  two 
in  phase  piston  sources  in  a  baffle  with  a  closed  ended  tube  between  them.  The 
tube  is  tuned  to  a  one  quarter  wavelength  resonance  of  200  Hz.  The  outer  sources 
are  separated  by  36cm,  center  to  center. 

This  configuration  was  devised  to  be  compared  with  the  case  of  three  point 
sources  with  the  center  source  driven  out  of  phase  of  the  outer  sources,  which  was 
discussed  in  Section  2.5.1.  At  200Hz,  Figure  2.4,  the  active  intensity  is  pointing 
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Figure  6.1  The  measurement  configuration  of  two  baffled  pistons 
with  a  resonating  tube  tuned  to  a  one-quarter  wave¬ 
length  resonance  of  200Hz. 


into  the  center  source,  indicating  that  it  is  absorbing  power.  However,  the 
reactive  intensity  is  pointing  out  of  a1!  three  sources  because  there  is  a  pressure 
maximum  above  each  point  of  vibration.  Consequently,  to  locate  these  three 
sources  the  reactive  intensity  must  be  used.  The  source  localization  capabilities 
of  the  reactive  intensity  are  shown  very  clearly  by  Degeorges  and  Tichy  (1986). 

The  two  source  and  tuned  resonator  configuration  was  intended  to  mirror  the 
situation  with  the  three  sources  when  the  center  source  appeared  to  be  absorbing 
power.  In  this  case,  we  have  a  resonator  playing  the  role  of  a  power  absorber. 
Intuitively  we  can  say  that  at  the  resonant  frequency  of  the  tube,  the  impedance 
of  the  tube  will  be  very  low  so  the  tube  will  absorb  energy.  However,  at  an  anti¬ 
resonant  frequency  the  impedance  will  be  very  high  at  the  tube  opening,  so  that 
the  energy  will  propagate  to  the  farfield.  Consequently,  the  analysis  is  confined 
to  the  first  two  resonant  frequencies,  200Hz  and  600Hz,  and  to  an  anti-resonance, 
400Hz. 

The  measurements  were  made  over  two  planes  for  each  frequency.  Measure¬ 
ments  were  conducted  in  a  vertical  plane  along  the  symmetric  cross  section  of 
the  sources  and  perpendicular  to  the  baffle.  Since  this  plane  is  symmetric  with 
respect  to  the  sources  only  the  two  components  of  the  intensity  within  the  plane 
needed  to  be  measured.  Additional  measurements  were  conducted  in  a  second 
plane  parallel  to  the  baffle  and  3cm  above  the  sources.  Measurements  were  made 
in  the  vertical  plane  to  investigate  the  radiation  to  the  farfield  and  the  source 
interaction,  and  in  the  horizontal  plane  to  investigate  the  characteristics  of  the 
sources  and  the  tube  and  their  boundaries  with  the  baffle. 

A  quick  overview  of  the  measured  data  reveals  that  the  data  are  acceptably 
symmetric,  as  would  be  expected  with  the  configuration.  Looking  at  the  wave 
fronts  we  can  see  that  *he  right  source  lags  the  left  source  by  about  3  degrees.  The 
same  signal  was  input  into  each  piston  source,  but  it  is  within  reason  for  them 


not  to  respond  identically.  Three  degrees  is  1%  of  a  cycle  which  is  acceptable. 
This  minor  delay  accounts  for  the  non  symmetry  in  the  active  intensity  data.  A 
similar  bias  is  seen  at  all  frequencies. 

6.2.1.  Fundamental  Resonance  of  the  Tube 

Let  us  first  look  in  detail  at  the  measurements  at  200Hz,  Figures  6.2  to  6.6. 
The  active  intensity  in  Figure  6.2  is  pointing  out  of  the  two  sources  and  into 
the  resonator.  The  resonator  is  absorbing  power  because  it  effectively  dissipates 
energy.  This  energy  comes  from  the  two  sources.  At  200Hz  the  sources  are 
strongly  coupled,  so  a  large  part  of  the  energy  is  consumed  by  the  resonant  tube. 
The  wave  front  propagation  shows  that  the  waves  combine  from  the  two  sources 
and  at  the  point  marked  by  a  peak  in  the  phase  speed,  Figure  6.3  ,  is  split  into 
waves  that  enter  the  tube  and  waves  which  propagate  to  the  farfield. 

The  reactive  intensity  plot  shows  some  fascinating  results.  These  are  sensible 
results,  and  demonstrate  the  powerful  analysis  capabilities  of  the  intensity 
analysis.  The  reactive  intensity  points  out  of  the  sources  since  there  is  a  pressure 
maximum  over  each  source.  However,  the  reactive  intensity  points  into  the 
resonating  tube.  Consequently  all  the  energy,  even  that  associated  with  the 
reactive  intensity  is  being  input  to  the  tube  to  drive  it  at  resonance. 

These  results  provide  a  good  comparison  between  power  absorbed  by  a 
passive  source  and  an  active  source.  In  Figure  2.6  the  active  intensity  is  directed 
into  the  center  source,  which  is  driven  out  of  phase  of  the  two  outer  sources. 
However,  the  reactive  intensity  in  Figure  2.6  points  out  of  the  active  center 
source  because  the  pressure  is  a  local  maximum  directly  over  the  source,  even 
though  there  is  a  time  averaged  net  flux  of  energy  into  the  center  source.  In  the 
case  of  the  resonating  tube,  Figure  6.2,  the  pressure  at  the  open  end  of  the  tube 
is  a  local  minimum  because  of  the  boundary  conditions  at  resonance,  so  that  the 
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Figure  6.2  The  measured  active  intensity  and  wave  fronts  and 
reactive  intensity  at  200Hz  in  a  plane  perpendicular 
to  the  baffled  containing  two  in  phase  piston  sources 
and  a  tube  with  a  200Hz  fundamental  resonance. 
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Figure  6.6  The  numerical  values  of  the  real  and  imaginary  parts 
of  the  impedance  normal  to  the  baffle  at  200Hz.  The 
measurement  plane  is  parallel  to  and  3cm  above  the 
baffle  containing  two  in  phase  piston  sources  and  a 
tube  with  a  200Hz  fundamental  resonance. 
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reactive  intensity  points  toward  the  tube  opening.  These  two  examples  show  that 
when  the  active  intensity  measurements  indicate  an  absorbtion  of  power,  as  in 
Figures  2.6  and  6.2,  the  direction  of  the  reactive  intensity  distinguishes  between 
an  active  source  coupled  with  other  active  sources,  Figure  2.6,  or  a  passive  source 
absorbing  energy  to  be  driven,  Figure  6.2. 

For  the  resonance  to  occur  in  the  tube  the  pressure  is  a  maximum  at  the 
closed  end  of  the  tube  and  zero  at  the  open  end  of  the  tube.  There  are  two 
effects  which  should  be  emphasized.  First  the  pressure  is  a  minimum  over  the 
tube  opening.  The  pressure  is  finite,  not  zero,  because  there  are  losses  in  the  tube 
and  the  boundary  conditions  at  the  tube  opening  are  not  exactly  pressure  release, 
consequently,  there  is  active  intensity  pointing  into  the  tube.  More  interestingly, 
note  that  the  pressure  minimum  forms  an  arc  over  the  tube  opening.  The  arc  is 
best  visualized  by  the  reactive  intensity.  This  pressure  minimum  is  not  directly 
at  the  tube  opening,  because  of  the  end  effects  at  the  tube  opening.  Fundamental 
theories  show  that  there  is  a  mass  of  fluid  attached  with  the  air  vibrating  in  the 
tube.  This  extends  the  effective  length  of  the  tube  beyond  its  physical  length. 
Consequently,  the  location  of  the  pressure  minimum  marks  the  effective  end  of 
the  tube. 

It  is  interesting  to  note  that  the  added  mass  is  not  a  cylindrical  extension  of 
the  tube,  but  has  a  hemispherical  shape,  Figure  6.3.  Assuming  symmetry,  which 
the  horizontal  plane  measurements  show,  the  volume  of  fluid  in  this  hemisphere 
is  |7r(.05)3  =  261cm3.  This  corresponds  to  an  equivalent  cylinder  of  fluid  1.5cm 
high.  The  theoretical  end  correction  for  an  open  tube  is  0.84  times  the  tube  radius 
if  flanged  and  0.6  times  the  radius  if  unflanged,  which  for  this  tube  is  6.6cm  and 
4.8cm  respectively.  In  both  predictions  an  end  correction  much  greater  than  the 
measurements  are  implied.  Most  likely,  this  is  a  consequence  of  the  acoustic 
coupling  between  the  two  piston  sources  and  the  tube.  The  coupling  alters  the 
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boundary  conditions  at  the  tube  opening  from  the  ideal  pressure  release.  The 
sources  and  tube  are  acoustically  separated  by  O.lA  so  the  tube  opening  is  exposed 
to  a  large  direct  field  from  the  sources,  which  alters  the  impedance  of  the  tube 
opening  and  thereby  reducing  the  extension  of  the  tube  length  beyond  the  end 
of  the  tube. 

Further  support  for  the  added  oscillating  mass  is  seen  in  the  kinetic  energy 
plots,  Figure  6.4.  There  is  a  region  above  the  tube  where  the  kinetic  energy 
is  large.  However,  this  is  confined  to  a  very  close  distance  above  the  tube 
opening,  compared  to  the  region  above  sources  which  extends  much  further. 
The  maximum  velocity  over  the  tube  extends  from  the  edge  of  the  tube  in  a 
semi-circle  to  about  5cm  above  the  tube  opening.  This  is  best  visualized  from 
the  plot  of  the  numeric  values. 

The  next  important  feature  of  this  measurement  is  the  normal  specific 
acoustic  impedance:  the  real  and  imaginary  parts,  Figures  6.5  and  6  6.  The 
measurements  made  on  the  plane  parallel  to  the  baffle  cover  the  area  over  the 
entire  surface  of  the  piston  sources  and  the  tube  opening.  The  real  part  of 
the  complex  normal  impedance  has  significant  minima  over  the  piston  sources, 
corresponding  to  power  radiation.  There  is  also  a  minimum  above  the  tube 
because  energy  is  being  absorbed.  Now  looking  at  the  imaginary  part  there  are 
weak  minima  over  the  two  piston  sources  because  there  is  a  reactive  component 
associated  with  their  radiation.  However,  there  is  a  deep  minimum  directly  over 
the  tube.  This  is  expected  above  a  resonator  and  is  graphically  shown  here.  In  a 
later  section  these  results  will  be  further  analyzed  by  comparing  them  with  the 
results  at  other  frequencies. 

A  curious  feature  of  the  measurement  is  the  large  active  intensity  pointing 
into  the  tube  opening,  Figure  6.2.  The  kinetic  energies  associated  with  the 
particle  velocities  Ua  and  Ur  are  roughly  of  the  same  amplitude  in  this  region. 
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Since  the  tube  is  a  resonator  one  might  expect  to  have  a  pure  standing  wave  in 
the  tube  so  that  the  kinetic  energy  is  all  contained  in  the  component  due  to  Ur ■ 
Consequently,  it  is  expected  that  the  kinetic  energy  just  above  the  tube  would 
be  dominated  by  the  component  due  to  Ur.  However,  a  more  realistic  analysis  of 
the  field  in  the  tube  must  include  damping  from  vibration  of  the  tube  walls  and 
air  gaps.  The  tube  which  is  being  driven  at  a  resonance  is  an  efficient  hoarder  of 
energy  so  there  is  a  large  response  in  the  tube.  Consequently,  any  small  damping 
in  the  tube  effectively  absorbs  energy  because  the  tube  easily  accepts  energy  and 
is  an  efficient  absorber  resulting  in  a  large  power  flux  into  the  tube. 

Another  important  aspect  of  the  impedance  above  the  tube  shown  in  Figures 
6.5  and  6.6  is  that  the  real  part  does  not  go  to  zero  because  of  the  energy  being 
absorbed  by  the  tube.  The  fundamental  equations  which  give  the  end  correction 
of  0.84  times  the  tube  radius  also  specify  that  the  real  part  of  the  impedance  at 
the  tube  opening  is  zero.  Consequently,  the  losses  in  the  tube  and  the  coupling 
with  the  sources  has  shifted  the  predicted  resonance  of  the  tube  with  the  extended 
length  by  decreasing  the  effective  length  of  the  tube.  This  is  done  by  shifting 
kinetic  energy  associated  with  the  added  mass  which  would  be  Ur ,  into  kinetic 
energy  due  to  the  component  Ua  which  inputs  power  into  the  tube  where  it  is 
dissipated. 

6.2.2.  Anti-Resonance  of  the  Tube 

Now  let  us  compare  the  results  at  resonance  with  those  at  the  first  anti¬ 
resonance  of  the  tube,  400Hz.  At  anti-resonance  the  closed  tube  has  a  very  high 
input  impedance,  so  that  it  is  essentially  not  seen  by  the  acoustic  field.  Therefore, 
the  field  should  be  much  the  same  as  that  for  two  baffied  pistons. 

The  active  intensity,  Figure  6.7,  clearly  shows  that  there  is  no  power  input 
to  the  tube  and  that  the  waves  propagate  over  the  tube  unaltered.  The  reactive 
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Figure  6.7  The  measured  active  intensity  and  wave  fronts  and 
reactive  intensity  at  400Hz  in  a  plane  perpendicular 
to  the  baffled  containing  two  in  phase  piston  sources 
and  a  tube  with  a  200Hz  fundamental  resonance. 


intensity  around  the  tube  opening  is  different  than  what  would  be  expected  than 
for  two  baffled  pistons  alone.  The  active  and  reactive  intensity,  with  the  tube 
having  an  infinite  input  impedance,  should  not  be  effected  by  the  tube.  The 
maximum  level  of  the  reactive  intensity  has  decreased  by  5dB  compared  to  the 
maximum  level  of  the  reactive  intensity  at  resonance  (200Hz)  This  decrease  is 
due  to  a  smaller  influence  of  the  tube  at  anti-resonance  and  the  lower  mutual 
coupling  between  the  pistons  at  this  higher  frequency. 

The  component  of  the  kinetic  energy  normal  to  the  baffle  is  plotted  in  Figure 
6.8.  The  kinetic  energy  normal  to  the  tube  opening  is  very  low.  Thus  the  tube 
is  acoustically  a  rigid  boundary,  as  indicated  in  the  previous  paragraph. 

The  normal  impedance,  Figure  6.9  and  6.10,  measured  in  the  plane  parallel  to 
the  baffle  shows  a  sharp  contrast  between  the  anti-resonance  and  first  resonance 
of  the  tube.  The  real  part  of  the  normal  impedance  at  400Hz  has  a  significant 
dip  at  the  sources.  However,  above  the  tube  the  real  impedance  has  the  same 
value  as  above  the  rigid  baffle.  The  one  peak  is  due  to  a  measurement  anomaly. 
These  results  show  that  the  power  radiated  at  this  frequency  is  only  a  function 
of  the  coupling  between  the  two  sources  and  that  the  tube  is  the  same  as  the 
rigid  baffle. 

The  imaginary  part  of  the  normal  impedance  has  slight  dips  over  the  piston 
sources,  which  is  associated  with  the  reactive  component  of  their  radiation. 
Above  the  center  of  the  tube  the  imaginary  part  is  the  same  level  as  over  the  baffle 
except  for  a  ring  where  there  is  a  local  minimum.  The  high  level  of  the  imaginary 
part  over  the  tube  agrees  with  the  high  level  of  the  real  part,  and  the  minimum 
ring  is  not  as  deep  as  the  imaginary  part  measured  at  resonance.  The  plot  of 
the  numeric  values  shows  that  the  dip  in  the  imaginary  part  occurs  at  the  edge 
of  the  tube.  Consequently,  we  can  conclude  that  there  is  some  slight  diffraction 
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occurring  because  of  the  sharp  edge  of  the  tube  opening.  This  diffraction  at  the 
tube  opening  causes  the  variation  in  the  reactive  intensity  seen  in  Figure  6.7. 

These  results  graphically  show  the  well  known  effects  of  using  resonators  to 
absorb  energy.  At  400Hz  the  damping  capabilities  of  the  tube  walls  is  no  different 
than  at  200Hz.  However,  at  400Hz  the  tube  has  a  very  high  input  impedance  so 
that  any  energy  in  the  vicinity  of  the  tube  would  rather  propagate  toward  the 
low  impedance  farfield  than  into  the  tube.  Without  any  energy,  the  damping 
within  the  tube  has  nothing  to  dissipate.  Therefore,  an  issue  in  absorbing  energy 
is  to  first  induce  the  energy  to  regions  with  damping. 

6.2.3.  Second  Resonance  of  the  Tube 

The  second  resonance  of  the  tube  occurs  at  600Hz.  The  field  is  much  the 
same  as  at  200Hz,  but  some  changes  are  very  important  to  notice,  Figure  6.11. 
Since  the  tube  is  operating  at  resonance  the  energy  is  once  again  being  absorbed 
by  the  tube.  At  resonance  there  is  sufficient  energy  in  the  tube  for  the  damping 
to  be  effective.  However,  the  distance  where  the  energy  being  absorbed  does  not 
extend  as  far  from  the  tube  as  at  200Hz. 

Despite  the  small  extent  from  the  tube  where  the  energy  points  in,  the  arc 
of  pressure  minima,  seen  clearly  in  the  reactive  intensity  in  Figure  6.11,  has 
the  same  size  as  at  200Hz.  This  is  expected  because  the  tube  length  correction 
at  the  lower  modes  is  independent  of  frequency.  It  is  especially  interesting  to 
compare  the  normal  kinetic  energy  at  600Hz,  Figure  6.12,  with  200Hz,  Figure 
6.4.  The  maximum  amplitude  is  lower  at  600Hz,  but  the  relative  distributions  of 
the  kinetic  energy  on  the  two  lines  of  data  above  the  baffle  are  the  same  at  these 
two  frequencies.  This  signifies  that  the  pistons  and  tube  are  constant  velocity 
sources,  which  are  not  altered  by  the  mutual  coupling  of  the  sources  and  the 
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The  specific  impedance  shown  in  the  plane  parallel  to  the  baffle,  Figure  6.13, 
has  the  same  characteristics  at  600Hz  as  at  200Hz.  The  dips  in  the  real  part 
are  confined  much  closer  to  the  sources  and  tube  opening  because  of  the  higher 
frequency,  which  results  in  the  lower  mutual  coupling.  In  fact,  this  is  apparent 
to  a  greater  extent  in  the  imaginary  part.  It  is  perhaps  easier  to  compare  the 
impedance  at  the  two  frequencies  by  using  the  numeric  values,  Figure  6.14.  Many 
of  the  features  of  these  plots  are  decreased  geometrically  in  size  by  about  3  times, 
corresponding  to  the  change  in  wavelength. 

6.2.4.  Comparison  of  the  Impedance  Above  the  Tube 

The  impedance  measured  above  the  tube  is  summarized  in  Table  6.1. 
Theoretically  at  the  tube  opening  the  impedance  is  purely  imaginary  and 
increases  linearly  with  frequency  at  resonance,  which  would  mean  a  4.8  dB  rise 
between  200  and  600  Hz.  The  imaginary  part  of  the  impedance  follows  this 
reasonably,  increasing  by  4  dB.  The  same  theory  predicts  that  the  real  part  will 
increase  with  the  square  of  the  frequency  or  9.5  dB.  However,  the  measured 
change  is  negligible.  One  explanation  of  this  is  that  at  resonance  the  damping 
in  the  tube  is  dominating  the  power  input  to  the  tube,  therefore  the  dominant 
influence  on  the  impedance.  Thus,  the  tube  accepts  energy  much  more  readily 
than  the  resonant  mechanism  would  predict.  The  other  explanation  is  mutual 
coupling,  which  is  frequency  dependent. 

6.3.  Variations  in  Radiation  Impedance  and  Output  Power 

The  measurements  which  are  reported  in  this  section  were  performed  at  the 
Physikalisch-Technische  Bundesanstalt  (PTB)  in  Braunschweig,  West  Germany. 
The  instrumentation  was  different  from  that  which  was  used  for  all  the  other 
measurements  in  this  thesis,  so  the  description  of  the  measurement  setup  is 
presented  in  Section  6.3.1. 
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Figure  6.14  The  numerical  values  of  the  real  and  imaginary  parts 
of  the  impedance  normal  to  the  baffle  at  600Hz  (the 
second  resonance  of  the  tube).  The  measurement 
plane  is  parallel  to  and  3cm  above  the  baffle  containing 
two  in  phase  piston  sources  and  a  tube  with  a  200Hz 
fundamental  resonance. 
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Table  6.1  The  real  and  imaginary  parts  of  the  impedance  normal 
to  the  baffled,  averaged  over  the  tube  opening.  The 
frequencies  200Hz  and  600Hz  are  at  the  tube  reso¬ 
nance  and  400Hz  is  at  the  tube  anti-resonance. 


Frequency 

Re{z } 

Im{z} 

(Hz) 

(dB) 

(dB) 

200 

13 

11 

400 

28 

16 

600 

14 

15 

The  output  power  of  a  simple  source  can  be  calculated  from  the  product  of 
the  source  volume  velocity  and  real  part  of  the  specific  acoustic  impedance  at  the 
surface  of  the  source.  Assuming  that  a  source  has  a  constant  volume  velocity, 
any  change  in  the  output  power  is  equated  to  the  variation  of  the  radiation 
impedance.  Variations  of  the  output  power  and  radiation  impedance  occur  when 
the  distance  of  a  source  above  a  reflecting  plane  changes.  Such  variations  are 
confined  to  low  frequencies  and  are  of  concern  when  determining  the  output 
power  of  a  source.  The  height  of  the  source  during  measurements  may  produce 
power  estimates  much  different  than  the  power  which  the  source  radiates  in  its 
normal  environment.  Therefore,  the  purpose  of  these  measurements  was  to  show 
that  the  intensity  technique  can  measure  the  impedance  at  the  surface  of  a  source 
which  can  then  be  compared  with  the  measurement  of  the  output  power. 

An  important  result  of  these  measurements  was  the  observation  that  signal 
processing  errors  in  the  real-time  analyzer  were  introduced.  The  error  was  due 
to  a  phase  shift  between  the  separate  analog  circuits  which  calculate  the  pressure 
and  velocity.  Thus,  errors  introduced  by  the  two  microphone  intensity  technique 
are  not  confined  to  the  phase  mismatch  between  the  microphones  and  the  two 


channels  of  the  signal  processing.  This  is  an  important  issue  for  institutes  which 
are  responsible  for  certifying  the  accuracy  of  equipment. 

There  is  some  notation  which  is  only  used  in  this  section.  The  notation  SPL 
will  be  used  for  the  sound  pressure  level,  SV L  for  the  sound  velocity  level,  SIL 
for  the  sound  intensity  level  (active  intensity),  and  SQL  for  the  reactive  intensity 
level.  This  is  done  to  comply  with  the  notation  used  at  the  PTB. 

6.3.1.  Measurement  Setup 

For  these  measurements  a  B&K  j  inch  phase  matched,  face-to-face  intensity 
probe  was  used.  The  small  size  was  needed  to  get  close  to  the  source  surface  while 
minimizing  diffraction  by  the  probe.  For  these  two  reasons  larger  devices  such  as 
the  B&K  ^  inch  probe  and  the  Norwegian  Electronics  intensity  probe  (Bjor  1985) 
were  not  used.  The  Norwegian  Electronics  intensity  probe  measures  the  pressure 
with  a  single  microphone  and  the  velocity  with  four  ultrasonic  transducers. 
The  ultrasonic  transducers  measure  doppler  shifts  and  are  orientated  on  a  line 
perpendicular  to  the  axis  of  the  microphone.  In  this  application  the  Norwegian 
Electronics  probe  also  has  the  disadvantage  of  being  sensitive  to  air  currents 
since  the  velocity  is  measured  directly.  Consequently,  with  a  wind  screen  and 
the  orientation  of  the  ultrasonic  transducers  and  the  microphone,  the  probe  was 
to  bulky  to  be  effective,  for  these  measurements. 

A  Norwegian  Electronics  one-third  octave  real-time  analyzer  type  830  was 
used  for  the  signal  processing.  This  is  a  two-channel  analyzer  which  measures 
in  one-third  octave  bands  the  sound  pressure  level  (SPL),  sound  velocity  level 
( SV L ),  and  sound  intensity  level  (SIL).  The  reactive  intensity  was  not  measured 
directly  by  the  analyzer,  so  it  needed  to  be  calculated  after  the  measurement 
was  completed.  The  active  intensity  (SIL)  and  reactive  intensity  (SQL)  are 


calculated  from  the  amplitudes  of  the  pressure  and  velocity  and  the  phase  angle 
6  between  them: 

I  =  cos6  ,  (6.1) 

Q  =  sin6  .  (6.2) 

2 

From  Equation  6.1  the  phase  angle  6  can  be  calculated: 

10  logio(cosO)  =  SIL-  -(SPL  +  SVL)  .  (6.3) 

mt 

Note  that  the  ^  in  Equation  6.1  is  not  included  in  Equation  6.3  because  SPL 
and  SV  L  are  effective  levels.  The  reactive  intensity  level  can  now  be  calculated: 

SQL  =  -( SPL  +  SVL )  +  10  logio(stnd)  ,  (6.4) 

A 

which  allows  the  real  and  imaginary  parts  of  the  impedance  to  be  calculated: 

10  /ogl0(/2e{z}}  =  SIL  -  SVL  ,  (6.5) 

10  log10(Im{z})  =  SQL  —  SV L  .  (6.6) 

The  measurements  were  performed  in  a  large  hemi-anechoic  room.  The 
floor  was  a  smooth  rigid  surface  and  the  walls  and  ceiling  were  anechoic.  For 
the  impedance  measurements  the  probe  was  positioned  by  hand.  The  power 
measurements  were  made  in  compliance  with  the  international  standards  for 
measuring  power  using  pressure  measurements;  however,  the  intensity  probe  was 
used.  The  probe  was  swept  in  a  quarter  of  an  arc  of  2m  in  radius.  The  intensity 
was  continuously  averaged  over  the  12^  minutes  of  this  process.  The  source  is 
rotated  so  the  path  of  the  probe  relative  to  the  source  traverses  a  hemisphere 


along  a  spiral  path.  The  source  which  was  used  was  omni-directional  so  that  the 
rotation  gave  similar  results. 

A  B&K  reference  sound  source  type  4205  was  used  for  the  measurement, 
Figure  6.15.  This  broad  band  noise  source,  which  uses  a  loudspeaker,  was  chosen 
over  the  fan  sources  because  of  the  large  air  motion  around  the  fan.  Also,  the  the 
radiation  surface  of  the  loudspeaker  is  small,  so  the  source  is  omni-directional 
and  can  be  classified  as  simple.  Therefore,  the  impedance  measured  at  one  point 
is  sufficient  to  characterize  the  source  impedance. 

The  acoustic  impedance  measurements  were  made  at  the  point  marked  on 
the  sketch  of  the  source,  Figure  6.15.  This  point  was  chosen  as  close  to  the  center 
of  the  radiating  surface  as  the  probe  could  be  positioned.  The  impedance  and 
radiated  power  was  measured  as  a  function  of  source  height  above  the  floor.  The 
source  was  raised  20cm,  40cm,  and  80cm  above  the  floor. 

6.3.2.  Analysis  of  the  Results 

The  real  and  imaginary  parts  of  the  specific  acoustic  impedance  measured  in 
one-third  octave  bands  centered  from  63Hz  to  5000Hz  is  shown  in  Figure  6.16. 
The  results  are  very  similar  to  those  calculated  for  a  small  sphere  with  ka  small, 
where  a  is  the  radius.  At  frequencies  below  800Hz,  where  ka  =  1.  the  impedance 
increases  with  frequency  and  the  imaginary  part  is  larger  than  the  real  part.  The 
real  and  imaginary  parts  are  equal  at  800Hz.  At  frequencies  above  800Hz  the 
real  part  of  the  impedance  stays  constant  and  the  imaginary  part  decreases.  In 
Figure  6.16  the  real  and  imaginary  parts  intersect  at  the  band  around  800Hz, 
which  corresponds  to  a  value  for  a  ~  7 cm.  This  is  only  1cm  smaller  than  the 
physical  size  of  the  source  radiating  surface.  These  results  show  that  the  source 
can  be  viewed  as  a  simple  source  so  that  the  impedance  can  be  measured  at  one 
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Figure  6.15  A  sketch  of  the  B&K  reference  sound  source  and  the 
position  of  the  measurement  point. 
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Next  the  source  was  raised  to  different  heights.  Figure  6.17  shows  the 
impedance  normalized  to  the  impedance  of  the  source  located  on  the  floor. 
Similarly,  the  radiated  power  as  a  function  of  the  source  elevation  referenced 
to  the  power  of  the  source  on  the  floor  is  plotted  in  Figure  6.18. 

The  overall  shape  of  these  curves  shows  the  change  in  the  radiated  power 
corresponding  to  the  change  in  the  radiation  impedance.  Above  250Hz  there  is 
excellent  agreement.  Below  250Hz  the  deep  minima  are  seen  in  both  the  power 
and  impedance,  but  the  change  in  impedance  is  much  less  than  the  corresponding 
changes  in  the  output  power. 

The  measured  SV L  varied  only  by  ±0.2 dB  for  each  source  height.  This 
variability  is  within  the  repeatability  of  the  measurement,  and  therefore  the 
source  can  be  considered  a  constant  velocity  source.  We  can  also  conclude 
that  the  measured  particle  velocity  amplitude  was  accurate  and  the  discrepancy 
between  the  power  and  impedance  changes  cannot  be  explained  by  the  variation 
in  the  SVL. 

The  discrepancy  could  also  have  been  caused  by  positioning  the  probe  to 
far  from  the  source  surface  so  that  the  measured  impedance  does  not  represent 
a  good  estimate  of  the  impedance  at  the  source  surface.  If  the  probe  was  not 
close  enough  then  we  would  expect  the  largest  discrepancies  at  high  frequencies 
where  the  distance  is  acoustically  large.  Likewise,  we  would  expect  the  most 
accuracy  at  low  frequencies  where  the  distance  is  acoustically  small.  However,  the 
measurements  above  250Hz  were  consistent  between  the  power  and  impedance 
so  the  probe  must  have  been  close  enough  below  250Hz.  Thus  the  distance  of 
the  probe  from  the  source  was  also  not  the  cause  of  the  discrepancy. 

The  only  other  significant  error  is  in  the  signal  processing  which  calculates 
the  sound  intensity  level.  Differences  between  the  SPL  and  SIL  of  lOdB  to 
15dB  are  beyond  the  dynamic  range  of  the  system.  But,  the  SPL  and  SIL 
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Figure  6.17  The  difference  in  dB  between  the  real  part  of  the 
impedance  measured  with  the  source  raised  20cm, 
40cm,  and  80cm  above  the  floor  and  the  source  on 
the  floor. 
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Figure  6.18 


The  difference  in  dB  between  the  radiated  power 
measured  with  the  source  raised  20cm,  40cm,  and 
80cm  above  the  floor  and  the  source  on  the  floor. 
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around  the  source  were  within  3dB  of  each,  so  the  intensity  level  was  not  beyond 
the  dynamic  range  of  the  system.  In  general,  phase  mismatch  between  the  two 
microphone  channels  causes  the  greatest  errors  at  low  frequencies.  However, 
the  power  was  measured  accurately  in  the  far  field,  so  the  system  can  measure 
the  SIL  accurately  below  250Hz.  Thus  the  error  was  not  caused  by  the  phase 
mismatch  between  the  two  microphone  signals. 

The  power  measurements  which  were  performed  in  the  farfield  differed  from 
the  impedance  measurements  which  were  performed  in  the  nearfield,  in  that  the 
level  of  the  velocity  out  of  phase  with  the  pressure,  Ur,  was  was  much  larger  in 
the  nearfield.  This  is  shown  in  Figure  6.16  where  below  250Hz  ’.he  imaginary 
part  of  the  impedance  is  5dB  to  lOdB  higher  than  the  real  part,  so  that  in  the 
nearfield  the  particle  velocity  is  nearly  ninety  degrees  phase  shifted  with  respect 
to  the  pressure.  This  corresponds  to  a  high  SQL.  In  contrast,  above  250Hz  the 
real  part  of  the  impedance  is  dominate  so  that  the  phase  between  the  pressure 
and  velocity  is  much  smaller  and  consequently  the  SIL  is  measured  accurately. 
Despite  the  high  level  of  the  SQL  in  the  nearfield  compared  to  the  farfield  the 
difference  between  the  SPL  and  SIL  remained  constant  in  the  nearfield  and  in 
the  farfield  and  was  only  ±1.5dB. 

Now  consider  how  the  analyzer  processes  the  signals,  Figure  6.19.  There 
are  two  separate  electrical  circuits  which  calculate  the  pressure  and  velocity. 
The  pressure  and  velocity  signals  are  subsequently  multiplied  to  form  the 
instantaneous  intensity  signal.  The  instantaneous  intensity  signal  is  then  time 
averaged  to  display  the  SIL.  However,  there  are  phase  shifts  introduced  between 
the  two  circuits  which  calculate  the  pressure  and  velocity.  This  phase  error  will 
have  different  effects  on  the  time  averaged  instantaneous  intensity  depending  on 
the  phase  between  the  pressure  and  velocity. 
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When  the  SQL  is  very  low  the  pressure  and  velocity  are  nearly  in  phase. 
Any  phase  error  between  the  circuits  calculating  the  pressure  and  velocity  will 
have  very  little  influence  on  the  SIL.  However,  when  the  SQL  is  very  high 
the  pressure  and  velocity  are  nearly  in  phase  quadrature,  so  small  phase  shifts 
introduced  by  the  two  circuits  will  greatly  influence  the  SIL.  This  is  exactly  the 
behavior  of  the  discrepancies  between  the  power  and  impedance  measurements. 

Tests  were  performed  with  the  analyzer  by  electrically  creating  signals  which 
simulated  high  SQL.  However,  these  were  done  only  with  a  phase  between 
pressure  and  velocity  of  fifty  degrees.  The  analyzer  accurately  displayed  the  SIL. 
However,  to  simulate  the  errors  encountered  during  the  impedance  measurements 
at  low  frequencies,  greater  phase  shifts  must  be  tested.  Thus,  when  making 
intensity  measurements  with  a  real-time  analyzer  concern  cannot  be  focused 
only  on  phase  mismatches  between  microphones.  The  electronic  circuits  which 
separately  calculate  the  pressure  and  velocity  must  also  be  accurately  phase 
matched  if  measurements  are  going  to  be  made  where  the  reactive  intensity  is 
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Chapter  7 


EXPERIMENTAL  ANALYSIS  OF  DIFFRACTION 


7.1.  Introduction 

The  details  of  diffraction  are  not  easily  understood  because  of  the  difficulty 
in  exploring  the  nearfield,  where  the  diffraction  occurs.  Acoustic  intensity  mea¬ 
surements  provide  a  detailed  visualization  of  diffraction  regions.  Consequently, 
the  jocation  of  diffraction  can  be  determined  and  the  effects  on  the  rest  of  the 
acoustic  field  can  be  seen.  Three  measurements  presented  in  this  chapter  will 
be  analyzed  in  a  qualitative  rather  than  quantitative  manner.  This  is  done  to 
emphasize  the  visualization  capabilities  of  the  active  and  reactive  intensity. 

7.2.  Diffraction  in  a  Loudspeaker  Nearfield 

Loudspeakers  are  commonly  analyzed  as  sources  of  sound.  However,  this 
chapter  will  focus  on  the  diffraction  in  the  nearfield  of  a  loudspeaker  because  it 
degrades  the  loudspeaker  radiation  characteristics.  The  desired  characteristic  of 
a  good  loudspeaker  is  a  uniform  acoustic  farfield.  However,  there  are  distortions 
of  the  sound  radiation  pattern  that  originate  in  the  nearfield,  thus  the  analysis 
should  be  performed  in  the  nearfield  where  the  sound  is  produced. 

The  loudspeaker  is  sketched  in  Figure  7.1.  It  was  designed  as  a  two  way 
system  with  a  6dB  per  octave  crossover  at  600Hz.  In  the  frequency  range  of  the 
crossover  the  acoustic  field  is  dominated  by  coupling  and  interference  of  the  sound 
radiated  by  the  woofer  and  tweeter.  The  measurements  which  will  be  discussed 
were  made  at  1500Hz  and  2500Hz.  This  is  far  above  the  crossover  frequency,  so 
only  the  tweeter  is  active.  Assuming  the  tweeter  radiates  like  a  piston  of  radius 
3cm,  then  ka  =  1.37  at  2500Hz.  Thus,  the  radiation  of  the  tweeter  is  essentially 
omni  directional  since  the  first  zero  in  the  main  lobe  occurs  at  ka  ^  3.5. 


Figure  7.1  A  sketch  of  the  loudspeaker  used  for  diffraction  mea 
surements. 


Two  measurement  planes  were  used.  A  plane  perpendicular  to  the  loud¬ 
speaker  surface  through  the  center  of  the  woofer  and  tweeter  was  chosen  to 
characterize  the  radiation  toward  the  farfield.  This  is  a  plane  symmetric  with 
respect  to  the  sources  so  only  two  dimensions  of  the  intensity  was  measured. 
Measurements  were  also  made  in  a  plane  parallel  to  the  loudspeaker  surface  to 
see  the  full  interaction  of  the  sources  and  box  edges.  This  required  a  three 
dimensional  measurement  of  the  intensity  which  was  made  3.5cm  above  the 
speaker  surface. 

Diffraction  will  intuitively  be  most  prevalent  at  impedance  discontinuities. 
The  most  obvious  areas  on  the  loudspeaker  are  the  outer  edges  of  the  box  where 
there  is  a  discontinuity  between  the  wood  surface  and  air.  The  other  significant 
area  is  the  edge  of  the  woofer.  The  woofer  edge  protrudes  1cm  above  the  box 
surface,  but  more  significantly  the  woofer  represents  an  abrupt  discontinuity 
between  the  thin  woofer  cone  and  the  rigid  wood  surface  of  the  box.  These 
measurements  identify  both  discontinuities. 

7.2.1.  Diffraction  Effects  at  1500Hz 

Let  us  begin  with  1500Hz.  The  active  intensity,  Figure  7.2,  shows  that  the 
wave  fronts  originate  from  the  tweeter  and  propagate  to  the  farfield  as  spherical 
waves.  There  is  some  distortion  of  the  waves  as  they  go  past  the  woofer.  The 
cause  of  the  distortion  is  easily  visualized  in  the  reactive  intensity. 

The  reactive  intensity  clearly  points  out  of  the  woofer  edge,  Figure  7.2; 
consequently,  there  is  a  pressure  maximum  at  the  woofer  edge,  Figure  7.3.  The 
kinetic  energy  associated  with  the  velocity  normal  to  the  speaker,  Figure  7.4,  also 
has  a  local  maximum  at  this  point.  The  reactive  intensity  in  the  horizontal  plane, 
Figure  7.5,  also  shows  the  pressure  maximum  above  the  woofer  edge.  Further, 
there  is  a  pressure  minimum  and  another  maximum  between  this  edge  of  the 
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Figure  7.3  The  potential  energy  measured  at  1500Hz  on  a  plane 
perpendicular  to  the  loudspeaker  surface. 
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Figure  7.4  The  kinetic  energy  due  to  the  particle  velocity  normal 
to  the  loudspeaker  surface,  measured  at  1500Hz  on  a 
plane  perpendicular  to  the  loudspeaker  surface. 
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woofer  and  the  tweeter.  Therefore  there  is  a  standing  wave  between  the  woofer 
edge  and  the  tweeter  caused  by  the  direct  wave  from  the  tweeter  interfering  with 
the  wave  reflected  by  the  woofer  edge. 


In  Figure  2.5  and  throughout  the  discussion  of  the  tube  and  piston  source 
measurements  in  Chapter  6,  we  saw  that  the  reactive  intensity  provided 
invaluable  insight.  Likewise,  in  these  measurements  the  reactive  intensity  in 
the  nearfield  reveals  more  important  details  than  does  the  active  intensity. 

The  effects  of  the  box  edge  can  also  be  seen  in  these  measurements.  This 
is  most  clearly  shown  in  the  measurements  on  the  plane  parallel  to  the  box, 
Figure  7.5.  The  active  intensity  is  not  significantly  effected  by  the  box  edge. 
However,  the  reactive  intensity  points  out  of  pressure  maxima,  3cm  inside  the 
edge,  between  the  tweeter  and  top  edge  and  the  tweeter  and  right  edge.  The 
normal  kinetic  energy  is  also  significant,  Figure  7.6.  There  are  regions  of  local 
velocity  maxima  just  inside  the  box  edge  which  are  due  to  the  diffraction  at 
the  box  edge.  These  areas  are  not  directly  above  the  edge  possibly  because  the 
measurement  plane  is  3.5cm  (0.15A)  from  the  box. 

The  normal  impedance  shows  some  especially  significant  results,  Figure  7.7. 
The  real  part  of  the  normal  impedance  is  clearly  reduced  by  3  to  4  dB  beyond 
the  box  edge.  This  is  a  consequence  of  the  change  between  the  high  impedance 
of  the  wood  and  the  low  impedance  of  the  air. 

The  diffraction  is  seen  as  a  5  to  7  dB  decrease  in  the  imaginary  part  of  the 
impedance  around  the  box  edge.  The  diffraction  does  not  create  or  sink  power, 
thus  there  is  little  change  in  the  real  impedance.  Consequently,  the  peaks  in  the 
kinetic  energy  near  the  box  edges  are  associated  with  the  reactive  intensity. 

It  should  also  be  noted  the  kinetic  energy  maxima  at  the  box  edges  are 
greatest  at  the  portion  of  the  edges  closest  to  the  tweeter,  Figure  7.6.  This  is 
also  seen  as  a  slightly  higher  imaginary  normal  impedance  at  the  same  points. 
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kinetic  energy  due  to  the  particle  velocity  normal 
he  loudspeaker  surface,  measured  at  1500Hz  on 
iane  parallel  to  and  3.5cm  above  the  loudspeaker 


This  will  be  more  effectively  illustrated  in  the  next  section  where  comparisons  of 
these  measurements  will  be  made  with  the  measurements  at  2500Hz. 

7.2.2.  Diffraction  Effects  at  2500Hz 

The  consequence  of  the  diffractions  at  1500Hz  do  not  significantly  alter  the 
wave  front  propagation,  but  is  confined  to  disrupting  the  pressure  distribution. 
At  2500Hz  the  same  diffraction  occurs,  but  has  a  significant  impact  on  the  wave 
front  propagation. 

First  consider  the  vertical  plane,  Figure  7.8.  The  wave  fronts  are  greatly 
distorted  in  a  region  extending  from  the  woofer  edge.  The  reactive  intensity  is 
once  again  pointing  out  of  the  woofer  edge  furthest  from  the  tweeter,  indicating 
a  pressure  maxima  due  to  the  diffraction  at  this  edge.  The  potential  energy, 
Figure  7.9,  has  a  ’eep  minimum  beginning  in  the  region  between  the  woofer  and 
the  tweeter  then  extending  to  the  farfield.  The  measurements  in  the  horizontal 
plane  clearly  show  the  diffraction  along  with  the  standing  wave  between  the 
woofer  edge  and  the  tweeter,  Figure  7.10.  However,  the  active  intensity  indicates 
that  there  is  a  wave  radiating  from  the  edge,  later  combining  with  the  wave  from 
the  tweeter.  The  wave  emanating  from  the  woofer  edge  interferes  with  the  waves 
from  the  tweeter  causing  the  band  of  low  pressure  seen  in  the  vertical  plane. 

One  explanation  for  the  region  of  low  pressure  could  be  the  beam  pattern  of 
the  tweeter.  However,  as  already  stated  a  piston  of  radius  3cm  would  have  the 
first  side  lobe  at  ka  =  3.5.  But  at  2500Hz  ka  =  1.37  so  no  side  lobes  in  the  beam 
pattern  can  occur.  Also  the  active  and  reactive  intensity  clearly  indicates  that 
the  interference  is  caused  by  the  diffraction  at  the  woofer  edge. 

The  kinetic  energy  over  the  surface  of  the  speaker  reveals  in  great  detail  the 
areas  of  diffraction,  Figure  7.11.  There  is  a  significant  peak  directly  above  the 
woofer  edge  and  there  are  peaks  at  the  box  edges  closest  to  the  tweeter.  At  this 


d I • tance  cm 


Maximum  L«v«l(dB):  70. 

i  a  a  000OOOOOOOOOOO0  o  o  a  • 

>  a  a  ooOOOOOOOOOOOOOOO  a  a 

>  o  o  a  000OOOOOOOOOOO00  o  a 

>  a  a  a  ooOOOOOOOOOOOOOo  o  a 

>oooo0OOO000O0OOO0oo  a  a 
»  o  o  a  a  0O000000OO0OOO  o  o  a 

>  a  o  o  ooOOOOOOOOOOOOOo  o  a  a 

■  aooooOOOOOOOOOOOOOoooo 

•  a  O  o  •  oOOOOOOOOOOOOO  O  O  O  O 
• -000000000000000000000 

*  • .oooOOOOOOOOOOOOOOOoo 

>a  . oooooOOOOOOOOOoooooo 
JO  -aoOOOO  OOQQOOO  O  O  O  O  O  O  O 

ooo.  oooooo  OOOOOOO  O  o  o  o  o  O  . 


)00  0  O  O  O  O  O  a 


39.3 

dl atone*  cm 


34.3 


Figure  7.9 


The  potential  energy  measured  at  2500Hz  on  a  plane 
perpendicular  to  the  loudspeaker  surface. 
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frequency  the  portions  of  the  edges  which  have  a  large  normal  velocity  are  shorter 
in  length  than  at  1500Hz.  This  seems  to  support  the  idea  that  the  greatest  effect 
of  the  edge  is  at  points  closest  to  the  tweeter  because  of  the  decrease  of  energy 
with  distance  from  the  tweeter.  Thus  the  diffraction  at  the  box  edge  does  not 
depend  strongly  on  the  incident  angle  with  the  edge. 

At  this  higher  frequency  the  impedance  is  less  usable  than  at  1500Hz,  Figure 
7.12.  This  is  partly  due  to  the  increased  variation  in  the  spatial  pressure 
and  velocity  distribution.  However,  at  2500Hz  the  plane  3.5cm  above  the 
speaker  is  0.26A  from  the  loudspeaker  surface,  as  opposed  to  0.15A  at  1500Hz. 
Consequently,  this  plane  is  acoustically  more  distant  from  the  surface,  thus 
imposing  a  limitation  for  the  nearfield  analysis.  If  the  specific  acoustic  impedance 
is  used  to  characterize  a  surface  then  the  measurement  must  be  acoustically  close. 
This  only  significantly  effects  the  impedance.  The  intensities  and  energy  densities 
still  provide  sufficient  information  at  this  distance  from  the  surface. 

7.3.  Nearfield  of  a  Circular  Disc 

This  example  presents  measurements  which  were  made  around  a  circular  disc, 
with  a  piston  source  1.5A  from  the  disc,  Figure  7.13.  Two  different  discs  were 
measured:  (l)  a  rigid  disc  constructed  of  |  inch  thick  particle  board  and  (2)  an 
absorptive  disc  constructed  of  8cm  thick  fiberglass.  The  measurements  were  made 
on  a  plane  perpendicular  to  the  disc  and  intersecting  the  disc  symmetrically.  The 
source  is  included  in  the  plane  so  that  the  intensity  does  not  have  components 
normal  to  the  measurement  plane.  A  more  complete  analysis  of  this  configuration 
would  include  a  three  dimensional  measurements  in  the  disc  nearfield,  especially 
over  the  surface  and  edge  of  the  disc.  However,  for  the  purposes  which  this 
example  serves  for  this  thesis,  the  analysis  is  confined  to  the  two-dimensional 
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Figure  7.13  Sketches  of  the  configurations  for  the  disc  diffraction 
measurements.  The  measurements  were  made  at 
1000Hz  on  a  plane  perpendicular  to  the  disc  and 
intersecting  the  source  and  disc  symmetrically. 


symmetric  plane.  The  frequency  is  1000Hz  so  that  the  relative  disc  size  is  1.22A 
and  the  measurement  extends  0.8A  from  the  disc. 

7.3.1.  Rigid  Disc 

This  first  measurement  demonstrates  the  beauty  of  the  acoustic  intensity 
measurement.  For  this  example  the  source  was  positioned  on  the  disc  axis 
so  that  a  full  rotational  symmetry  was  achieved.  The  active  and  reactive 
intensity  are  plotted  in  Figure  7.14  and  the  potential  energy  and  phase  speed 
are  plotted  in  Figure  7.15.  First  look  at  the  active  and  reactive  intensity.  There 
is  no  analytical  solution  to  this  source  scatterer  configuration  in  the  nearfield. 
However,  the  active  intensity  shows  the  wave  front  propagation  very  clearly.  The 
waves  approach  the  disc  then  divert  around  the  disc  circling  two  points  where 
the  pressure  is  zero.  These  two  formations  are  vortices  which  were  explained 
in  Section  4.5.  The  waves  then  combine  on  the  shadow  zone  of  the  disc  and 
propagate  nearly  as  plane  waves  to  the  farfield. 

The  reactive  intensity  indicates  a  standing  wave  between  the  source  and  disc 
which  arises  from  interference  of  the  direct  wave  of  the  source  and  the  wave 
reflected  by  the  surface  of  the  rigid  disc.  There  is  also  a  significant  pressure 
gradient  at  the  upper  and  lower  edge  of  the  disc  facing  the  source.  In  the  shadow 
zone  there  is  a  standing  wave  which  consists  of  a  horizontal  pressure  maximum 
and  minimum  caused  by  the  interference  of  energy  coming  over  the  top  and 
bottom  of  the  disc.  Consequently,  the  reactive  intensity  shows  that  the  region  on 
the  source  side  of  the  disc  has  a  vertical  standing  wave  and  in  the  shadow  zone 
has  a  radial  standing  wave. 

It  was  noted  that  there  is  a  pressure  maximum  at  the  corners  of  the  disc 
edge  facing  the  source.  The  pressure  maximum  is  seen  in  the  reactive  intensity 
which  points  away  from  these  corners.  The  area  to  focus  on  is  outlined  on  the 
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reactive  intensity  plot,  Figure  7.14.  The  edge  represents  a  sharp  discontinuity 
between  the  rigid  wood  and  the  air.  These  measurements  indicate  that  the  effect 
of  the  discontinuity  is  most  significant  on  the  corners  of  the  edge  facing  the 
source,  compared  to  the  corners  of  the  edge  facing  the  shadow  zone.  They  also 
indicate  that  the  wave  propagation  is  not  as  dramatically  effected  as  the  pressure 
distribution  so  that  the  reactive  intensity  most  readily  identifies  this  diffraction 
effect.  In  essence,  the  kinetic  energy  associated  with  the  diffraction  is  that  which 
is  due  to  the  component  of  velocity  which  is  out  of  phase  with  pressure,  Ur. 
To  control  this  diffraction  in  the  nearfield  it  seems  that  the  abrupt  impedance 
change  at  the  edge  facing  the  incident  sound  needs  to  be  reduced  with  absorption 
or  by  rounding  the  edge.  This  will  have  more  effect  than  treating  the  edge  in  the 
shadow  zone. 

Measurements  were  also  made  with  the  source  raised  to  the  level  of  the 
disc  top.  The  active  and  reactive  intensity  are  plotted  in  Figure  7.16  and  the 
potential  energy  and  phase  speed  are  plotted  in  Figure  7.17.  As  with  the  previous 
configuration  there  is  a  vertical  standing  wave  between  the  source  and  the  disc 
with  a  wave  circulating  a  pressure  zero  point.  Also,  in  the  shadow  zone  there  is 
a  horizontal  standing  wave  which  is  shifted  because  of  the  new  source  position. 

Once  again,  the  diffraction  at  the  corners  of  the  edges  facing  the  source  is 
clearly  visible.  The  diffraction  is  best  compared  with  the  case  of  the  symmetric 
source  by  using  the  acoustic  impedance.  Figure  7.18  contains  the  numerical 
values  of  the  real  and  imaginary  parts  of  the  impedance  close  to  the  disc  surface 
for  the  symmetric  and  raised  source  position.  The  real  part  of  the  impedance 
just  beyond  the  disc  edge  is  1  dB  lower  compared  to  the  level  at  the  disc  edge. 
For  the  symmetric  source  the  imaginary  part  of  the  impedance  just  beyond  the 
disc  edge  is  3  to  4  dB  lower  compared  to  the  level  at  the  disc  surface.  Similarly, 
with  the  raised  source  the  imaginary  part  of  the  impedance  decreases  by  4  to  6 
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Figure  7.18 


The  real  and  imaginary  parts  of  the  total  impedance 
around  the  disc.  The  top  figures,  Za,  are  for  the 
rigid  disc  with  the  symmetric  source  position  and  the 
bottom  figures,  Zr,  are  for  the  rigid  disc  with  the 
raised  source  position.  The  numbers  represent  the 
dB  level  below  the  maximum  level  labeled  above  each 
plot. 


dB  beyond  the  top  and  2  to  3  dB  beyond  the  bottom.  These  results  further  show 
that  the  imaginary  part  of  the  impedance  is  greatly  influenced  by  diffraction. 
Also  there  is  more  diffraction  at  the  top  of  the  disc  than  the  bottom  when  the 
source  is  raised  to  the  level  of  the  disc  top,  and  the  diffraction  with  the  symmetric 
source  is  between  these  two  extremes. 

The  distance  of  the  source  to  each  edge  is  one  possible  cause  of  the  variation 
of  the  diffraction  at  the  two  edges  for  the  two  source  positions.  When  the  source 
is  raised  to  the  level  of  the  top  edge,  the  distance  of  the  source  to  the  top  edge  is 
51.5cm  and  to  the  bottom  edge  is  65.2cm.  With  the  symmetric  source  position 
the  distance  to  the  two  edges  is  55.2cm.  The  amount  of  diffraction  increases 
as  the  source  to  edge  distance  decreases.  However,  this  is  not  a  likely  cause  of 
the  variation  in  diffraction.  The  various  source  distances  would  only  change  the 
amplitude  of  the  incident  wave;  however,  the  impedance  is  not  effected  by  the 
amplitude  of  the  incident  energy. 

The  appropriate  explanation  for  the  variation  in  the  diffraction  is  the  angle 
of  incidence  of  the  incident  wave,  which  is  the  other  major  geometric  variable 
in  these  measurements.  The  incident  angle  is  90  degrees  with  respect  to  the 
top  edge  with  the  source  raised,  69  degrees  with  respect  to  both  edges  with  the 
symmetric  source  position,  and  52  degrees  with  respect  to  the  bottom  edge  with 
the  source  raised.  Thus,  the  diffraction  increases  as  the  incident  angle  increases 
toward  90  degrees.  Physically  this  corresponds  to  a  greater  proportion  of  the 
incident  energy  being  normal  to  the  impedance  discontinuity  between  the  rigid 
disc  and  the  air.  These  results  are  identical  to  those  for  the  diffraction  at  the 
outer  edge  of  the  loudspeaker  box. 

This  explanation  can  be  further  verified  by  comparing  the  direction  of  the 
active  and  reactive  intensity  at  the  various  edges.  First  consider  the  top  edge 
for  the  raised  source.  The  wave  fronts  propagate  past  the  edge,  thus  the  active 


intensity  is  perpendicular  to  the  disc  plane  even  above  the  disc  top.  However,  the 
reactive  intensity  points  away  from  the  corner  of  the  top  edge  facing  the  disc,  thus 
the  energy  concentration  created  by  the  edge  discontinuity  must  propagate  away 
from  the  edge  via  the  reactive  intensity.  For  the  symmetric  source  configuration 
there  is  a  greater  portion  of  the  total  energy  parallel  to  the  disc  surface  contained 
in  the  active  intensity  rather  than  the  reactive  intensity.  Then  with  the  raised 
source,  the  active  intensity  at  the  bottom  edge  is  almost  parallel  to  the  disc 
surface  so  there  is  less  energy  associated  with  the  reactive  intensity  due  to  the 
diffraction. 

7.3.2.  Fiberglass  Disc  in  an  Anechoic  Room 

The  diffraction  by  a  fiberglass  disc  was  measured  to  show  a  strong  contrast  to 
the  rigid  disc  measurements.  The  fiberglass  disc  was  8cm  thick  and  of  the  same 
diameter  as  the  rigid  disc.  The  source  was  situated  along  a  line  extending  from 
the  top  of  the  disc.  The  other  parameters  such  as  the  frequency,  source  distance, 
and  measurement  plane  were  the  same  as  with  the  rigid  disc.  The  active  and 
reactive  intensity  are  plotted  in  Figure  7.19  and  the  potential  energy  and  the 
phase  speed  are  plotted  in  Figure  7.20. 

First  consider  the  region  between  the  disc  and  source.  The  reactive  intensity 
shows  a  standing  wave,  which  has  a  minimum  closer  to  the  disc  than  the  rigid 
case  because  of  the  impedance  of  the  fiberglass  disc.  There  are  no  pressure  zeros 
so  the  active  intensity  does  not  indicate  any  vortices,  but  shows  the  power  flux 
and  wave  propagation  into  the  disc.  Compared  with  the  rigid  disc,  the  active 
intensity  level  remains  approximately  the  same,  but  the  reactive  intensity  is 
reduced  by  3  dB.  It  is  interesting  to  note  that  the  partial  standing  wave  between 
the  source  and  the  disc  is  still  present  because  the  fiberglass  does  indeed  reflect 


some  energy. 
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Figure  7.20  The  potential  energy  and  resultant  phase  speed  mea¬ 
sured  at  1000Hz  on  a  plane  perpendicular  to  the  8cm 
thick  fiberglass  disc.  The  source  is  1.5A  to  the  left, 
along  a  line  extending  from  the  disc  top. 
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At  the  disc  edge  facing  the  source  there  is  a  diffraction  which  appears  to 


be  the  same  as  that  seen  with  the  rigid  disc.  However,  the  reactive  intensity 


indicates  that  the  diffraction  does  not  extend  as  far  from  the  edge  as  with  the 


rigid  disc.  This  is  expected  because  of  the  lower  impedance  of  the  fiberglass. 


The  numerical  values  of  the  real  and  imaginary  part  of  the  impedance  is  plotted 


in  Figure  7.21.  There  is  a  rapid  change  in  the  imaginary  part  of  the  impedance; 


therefore,  the  pressure  maxima  at  the  edge.  However,  this  rapid  change  is  now 


only  2  to  3  dB  which  is  smaller  than  the  4  to  6  dB  seen  with  the  rigid  disc. 


In  the  shadow  zone  the  data  is  greatly  different  for  the  fiberglass  and  rigid 


disc.  This  is  most  apparent  in  the  reactive  intensity  which  points  into  the 


fiberglass.  Consequently,  there  is  a  pressure  minimum  at  the  surface  of  the 


fiberglass,  and  no  horizontal  standing  wave.  This  is  due  to  the  attenuation  of 


the  energy  propagating  around  the  disc.  The  only  significant  energy  propagation 


parallel  to  the  disc  surface  is  near  the  edges  and  is  seen  only  in  the  reactive 


intensity.  The  active  intensity  is  pointing  out  of  the  disc,  which  indicates  that 


much  of  this  energy  comes  from  waves  propagating  through  the  disc. 


There  is  a  large  pressure  and  active  intensity  minimum  at  the  face  of  the 


disc  slightly  below  the  center  in  the  shadow  zone.  However,  there  is  a  significant 


amount  of  reactive  intensity  so  the  pressure  gradient  is  high.  This  pressure 


minimum  shows  that  the  energy  is  highly  attenuated  as  it  propagates  around  the 


disc.  The  minimum  is  below  the  disc  center  because  there  is  more  energy  coming 


over  the  top  of  the  disc  which  is  closer  to  the  source.  It  is  very  interesting  to 


note  that  at  this  region  where  the  pressure  and  active  intensity  have  a  low  level, 


the  mechanism  for  the  energy  removal  is  the  particle  velocity  associated  with  the 


the  reactive  intensity  which  points  into  the  fiberglass.  This  is  reasonable  if  we 


remember  that  there  is  energy  propagation  even  if  the  active  intensity  is  zero,  as 


was  shown  in  Chapter  2. 
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Figure  7.21  The  real  and  imaginary  part  of  the  total  impedance 
around  the  fiberglass  disc  with  the  raised  source 
position.  The  numbers  represent  the  dB  level  below 
the  maximum  level  labeled  above  each  plot. 


7.3.3.  Fiberglass  Disc  in  a  Reverberant  Field 

The  measurements  with  the  fiberglass  were  repeated  with  all  the  wedges 
removed  from  the  room,  thus  creating  a  reverberant  field.  The  reason  for 
performing  this  measurement  was  to  see  what  would  happen  with  the  energy 
associated  with  the  active  intensity  which  propagates  away  from  the  disc  under 
anechoic  conditions.  The  active  and  reactive  intensity  are  plotted  in  Figure  7.22 
and  the  potential  energy  and  resultant  phase  speed  are  plotted  in  Figure  7.23. 

Let  us  focus  on  two  aspects  of  this  measurement.  First  note  that  there  is 
no  significant  diffraction  at  the  disc  edge.  In  the  examples  of  the  rigid  disc  we 
saw  that  the  extent  of  diffraction  at  the  top  and  bottom  edge  depended  on  the 
incident  angle  of  the  energy.  In  the  reverberant  room  the  energy  is  incident  on 
the  fiberglass  disc  edges  at  all  angles  so  the  diffraction  is  reduced. 

The  other  important  aspect  of  this  measurement  is  the  direction  of  the  active 
intensity  which  points  into  the  fiberglass.  The  reactive  intensity  indicates  a 
build  up  of  local  energy  due  to  the  large  spatial  variation  of  the  pressure  in  the 
reverberant  conditions.  Consequently,  in  these  measurements  the  mechanism  for 
energy  propagation  to  the  fiberglass  resulting  in  energy  dissipation  is  the  active 
intensity. 

7.4.  Distinguishing  Between  Radiation  and  Diffraction 

This  last  example  investigates  the  acoustic  field  above  a  circular  plate  in 
the  shadow  zone  of  a  rigid  disc.  The  configuration  is  sketched  in  Figure  7.24. 
The  plate  is  constructed  of  ^  inch  thick  aluminum  and  is  fixed  at  its  edge.  The 
plate  was  intended  to  be  a  rigid  surface  so  that  the  diffraction  around  the  rigid 
disc  and  the  plate  and  supporting  structure  could  be  studied.  However,  the 
plate  can  be  excited  by  the  incident  sound.  The  ensuing  sound  radiated  by  the 
plate  vibration  would  disrupt  the  intended  of  investigation  of  diffraction  alone. 
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CIRCULAR  PLATE  36.5  cm  IN  DIAMETER 


Figure  7.24 


A  sketch  of  the  symmetric  cross  section  of  the  mea¬ 
surement  configuration  with  circular  plate  and  sup¬ 
porting  structure  in  the  shadow  zone  of  a  rigid  disc. 
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Consequently  the  following  measurements  were  made  to  indicate  the  influence  of 
the  plate  vibration  on  the  acoustic  nearfield  of  the  diffraction. 

The  plate  and  supporting  structure  was  placed  near  a  loudspeaker  radiating 
pink  noise.  The  acceleration  on  the  plate  and  the  supporting  structure,  measured 
with  an  accelerometer,  is  plotted  in  Figure  7.25.  The  resonant  response  of  the 
plate  is  20dB  to  30dB  larger  than  the  vibration  of  the  supporting  structure. 
Therefore,  attention  can  be  focused  on  the  resonances  of  the  plate  because  of 
the  weak  coupling  between  the  plate  and  supporting  structure.  Constrained 
layer  damping  was  applied  to  the  plate,  which  reduced  the  vibration  level  at 
resonance  frequencies  by  about  5dB. 

To  estimate  the  sound  radiated  by  the  vibration  of  the  plate  with  the 
constrained  layer  damping,  a  crude  measurement  was  made  with  the  plate 
positioned  behind  the  rigid  disc,  Figure  7.24.  By  placing  a  hand  on  the  plate 
the  vibration  levels  at  resonance  could  be  xeduced  by  10  dB.  The  pressure  2cm 
above  the  top  plate  was  measured  with  a  hand  touching  the  plate.  Then  the 
hand  was  lifted  off  the  plate  but  with  the  hand  in  essentially  the  same  location 
so  the  diffraction  was  not  altered.  The  small  variation  in  pressure  indicated  that 
the  radiated  pressure  at  a  resonance  of  the  plate  was  10  to  15  dB  below  the 
level  of  the  diffracted  field.  The  intensity  measurements  were  performed  at  anti¬ 
resonances  of  the  top  plate  vibration  so  it  was  considered  safe  to  assume  that 
the  sound  radiated  by  the  sphere  vibration  is  20  to  25  dB  below  the  level  of  the 
diffracted  sound  which  we  wanted  to  measure. 

However,  the  intensity  measurements  lead  us  to  believe  that  the  vibration  of 
the  plate  is  indeed  significant.  This  question  is  raised  by  the  results  at  1930Hz, 
Figure  7.26,  and  2570Hz,  Figure  7.27.  These  figures  show  an  expanded  section 
above  the  plate.  The  measurements  were  made  on  the  symmetric  axis  of  the 
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Figure  7.26 


An  expanded  view  of  the  measurement  over  the  plate 
in  the  shadow  zone  of  a  rigid  disc.  The  frequency  is 
1930Hz  and  the  source  is  53cm  (3A)  to  the  left  of  the 
disc  top. 
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source,  disc,  and  truncated  sphere  which  is  a  plane  perpendicular  to  the  face  of 
the  rigid  disc  and  includes  the  source  and  slices  the  truncated  sphere  in  half. 

Careful  inspection  of  the  reactive  intensity  just  above  the  top  plate  reveals 
areas  of  pressure  maxima  and  minimum.  These  regions  are  circled  in  the  plots. 
The  reactive  intensity  has  been  proven  useful  by  Degeorges  and  Tichy  (1986), 
to  locate  areas  of  maximum  vibration  on  rectangular  plates  so  that  the  mode 
shape  of  the  plate  is  readily  visualized.  In  general,  areas  of  maximum  pressure 
directly  above  a  vibrating  structure  will  correspond  to  areas  of  maximum  surface 
vibration.  Thus,  the  reactive  intensity  pattern  which  is  circled  can  be  interpreted 
to  indicate  the  modes  of  vibration  of  the  top  plate.  The  pattern  is  very  symmetric 
with  respect  to  the  top  plate  so  it  seems  very  possible  that  a  mode  of  the  plate 
is  significantly  excited. 

An  alternative  cause  for  this  pressure  distribution  is  to  consider  a  simple 
model  of  the  energy  propagating  over  the  top  of  the  disc  interfering  with  the 
energy  reflected  by  the  top  plate.  This  is  schematically  drawn  in  Figure  7.28. 
This  model  does  explain  the  pressure  minima  and  maxima  which  are  circled 
in  the  potential  energy  plots.  However,  the  paths  are  not  consistent  with 
half  wavelength  distances  which  would  cause  the  pressure  reinforcement  and 
cancellation  which  the  reactive  intensity  indicates  directly  above  the  top  plate. 
Thus  this  simplified  explanation  is  inadequate. 

The  definitive  measurement  which  should  have  been  made  was  on  a  plane 
parallel  to  the  top  plate  surface,  which  would  have  been  a  three  dimensional 
measurement.  This  should  indicate  a  mode  shape  of  the  top  plate  if  there  is 
significant  vibration.  However,  this  measurement  could  not  be  made  for  technical 
reasons.  If  the  measurement  indicates  a  mode  shape  then  the  active  intensity 
can  be  used  to  calculate  the  power  radiated  by  the  top  plate.  For  a  preliminary 
indication  Figures  7.2b  ~ud  7.27  can  be  used.  At  1930Hz  the  active  intensity 
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Figure  7.28  A  schematic  for  interference  of  a  direct  wave  over 
the  disc  top  and  a  wave  reflected  from  the  truncated 
sphere  top  plate. 


S.'Vjv 


145 


is  parallel  to  the  plate  so  there  is  only  a  negligible  power  flux  from  the  plate. 
Likewise  at  2570Hz  a  small  amount  of  power  goes  into  the  plate,  which  would 
indicate  that  the  plate  is  absorbing  rather  then  radiating  power.  However  at 
2570Hz  the  closest  measurement  point  is  ^  from  the  plate  surface  so  this  may 
not  be  a  valid  indication  of  the  power  flux  at  the  plate  surface.  These  results 
lead  to  the  conclusion  that  the  plate  is  radiating  no  power. 

Before  concluding  that  the  reactive  intensity  pattern  is  due  to  a  diffraction 
effect,  we  must  consider  the  nature  of  the  plate  vibration.  The  coincidence 
frequency  of  this  plate  is  roughly  19,000Hz.  The  measurement  frequencies  are  far 
below  the  coincidence  frequency,  so  the  plate  is  radiating  very  little  power.  Any 
acoustic  energy  radiated  by  the  plate  would  be  confined  to  the  nearfield  reactive 
intensity.  This  analysis  agrees  with  the  measurements.  The  radiated  pressure 
may  also  be  very  small  so  that  the  simple  measurement  described  earlier  in  the 
section  would  not  indicate  the  extent  of  the  plate  vibration.  However,  as  we  have 
seen  throughout  this  work,  the  reactive  is  sensitive  to  small  pressure  variations, 
so  that  in  these  measurements  the  reactive  intensity  indicates  the  extent  of  the 
non-radiating  nearfield  of  the  plate  vibrations.  Therefore,  by  considering  the 
reactive  intensity  we  can  conclude  that  the  plate  is  being  excited  by  the  incident 
sound  field  and  cannot  be  assumed  rigid. 
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Chapter  8 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
8.1.  Summary  and  Conclusions 

To  effectively  utilize  the  acoustic  intensity  in  the  nearfield,  the  physical 
interpretation  of  wave  phenomena  must  be  extended  beyond  the  simple  farfield 
or  plane  wave  properties.  The  conceptions  of  wave  motion  are  equally  valid  in 
the  nearfield,  but  the  wave  shape  and  propagation  speed  vary.  Consequently,  a 
description  of  the  energy  propagation  must  include  phenomena  other  than  the 
motion  of  the  wave  front. 

A  general  form  was  assumed  for  the  pressure  so  that  the  intensity  vectors 
could  be  linked  to  physical  phenomena.  The  inter-dependance  of  the  pressure 
amplitude  and  phase  distribution  are  expressed  by  two  non-linear  differential 
equations.  Since  the  waves,  in  general,  travel  at  speeds  other  than  c,  there  is  a 
non-uniform  pressure  amplitude  distribution  which  is  described  by  the  differential 
equations.  The  pressure  gradient  signifies  that  there  is  reactive  intensity,  thus  the 
wave  structure  moves  at  speeds  other  than  c,  because  there  is  energy  propagating 
associated  with  the  reactive  intensity. 

The  energy  propagation,  in  the  direction  of  the  instantaneous  intensity,  is 
in  general  not  perpendicular  to  the  wave  fronts.  Consequently,  in  general,  the 
energy  does  not  propagate  with  the  wave  fronts.  This  was  clearly  shown  with 
instantaneous  intensity  flux  lines  with  paths  dramatically  different  than  the  wave 
front  propagation.  However,  it  is  not  practical  to  use  the  instantaneous  intensity 
to  investigate  a  sound  field,  because  there  is  not  one  path  on  which  the  energy 
propagates  at  all  times. 
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The  time-independent  intensity  vectors  were  linked  to  fundamental  concepts 
of  wave  motion.  The  active  intensity  is  perpendicular  to  the  wave  fronts  and  the 
general  phase  speed  is  calculated  from  the  active  intensity  and  potential  energy. 
Contrary  to  the  farfield,  the  wave  shape  and  speed  changes  as  it  propagates.  The 
pressure  amplitude  variation  in  space  is  described  by  the  reactive  intensity  which 
points  in  the  direction  of  decreasing  pressure. 

The  active  intensity  vortex  is  an  interesting  example  which  combines  the 
concepts  of  energy  propagation  and  wave  front  propagation.  The  vortex  is 
centered  at  a  point  or  line  of  zero  pressure  which  also  corresponds  to  a 
discontinuity  in  the  pressure  phase.  At  the  phase  discontinuity,  different  surfaces 
of  constant  phase  intersect,  so  the  vortex  represents  a  wave  propagating  around 
the  line  or  point  of  zero  pressure.  Three  arguments  were  presented  which 
concluded  that  the  energy  propagates  through  the  vortex  rather  than  being 
trapped  by  the  vortex. 

The  capability  of  using  the  intensity  measurements  to  analyze  sound  fields 
is  the  most  important  reason  for  developing  the  physical  interpretations  of  the 
intensity  vectors.  The  five  measured  examples  explored  several  fundamental 
phenomena,  i.e.  sources,  absorption,  diffraction,  source-resonator  interactions, 
and  the  extended  length  of  a  resonating  tube.  The  fundamentals  of  the 
phenomena  are  well  understood.  In  the  past  these  phenomena  were  not  measured 
directly,  but  they  were  implied  to  be  the  cause  of  the  results  observed  in 
the  farfield.  However,  the  nearfield  measurements  in  this  thesis  gave  a  clear 
visualization  of  the  phenomena  and  showed  how  the  intensity  measurements  can 
be  used  to  investigate  more  complicated  phenomena. 

The  first  measurement  investigated  the  interaction  of  two  pistons  in  a  baffle 
with  a  resonating  tube  between  them.  The  tube  was  tuned  to  a  one-quarter 
wavelength  resonance  of  200Hz.  Measurements  were  performed  at  two  of  the 
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tube  resonances,  200Hz  and  600Hz,  and  at  a  tube  anti-resonance,  400Hz.  The 
significant  conclusions  from  these  measurements  were  (1)  the  reactive  intensity 
gave  a  visualization  of  the  extended  tube  length,  (2)  a  deep  minimum  in  the 
imaginary  part  of  the  normal  impedance  identified  the  tube  resonance,  (3)  sharp 
minima  in  the  imaginary  part  of  the  normal  impedance  identified  diffraction 
around  the  tube  edge  at  the  tube  anti-resonance,  and  (4)  a  through  analysis 
allowed  absorption  to  be  differentiated  from  sources. 

The  second  source  example  investigated  the  variation  in  the  radiated  power 
of  a  simple  source  as  a  function  of  the  height  of  the  source  above  a  rigid  boundary. 
The  measured  radiation  impedance  basically  agreed  with  the  change  in  the 
radiated  power.  The  few  discrepancies  which  were  observed  between  the  power 
and  impedance  measurements  were  due  to  the  pressure  and  velocity  being  nearly 
ninety  degrees  out  of  phase,  which  corresponds  to  a  high  reactive  intensity.  It 
turned  out  that  the  measurement  error  was  not  due  to  the  phase  errors  in  the 
probes,  but  due  to  the  phase  errors  between  the  processing  of  the  pressure  and 
velocity  in  the  real  time  one-third  octave  analyzer  which  was  used  for  this  one 
measurement. 

This  measurement  was  significant  in  two  aspects  because  it  pointed  out  the 
potential  problem  with  using  this  type  of  instrumentation.  A  real  time  analyzer 
will  often  only  display  the  pressure,  active  intensity,  and  the  particle  velocity. 
However,  the  reactive  intensity  and  the  real  and  imaginary  parts  of  the  impedance 
can  be  calculated  from  this  limited  information.  The  second  important  aspect 
was  that  with  a  real  time  analyzer  the  separate  electric  circuits  which  calculate 
the  pressure  and  velocity  need  to  be  carefully  phase  matched.  This  is  not  a 
problem  with  FFT  processing,  and  is  only  a  concern  with  real  time  analyzers  in 
the  presence  of  high  reactive  intensity  levels.  This  is  significant  because  intensity 
measurements  are  most  often  made  near  to  a  sound  source  where  the  reactive 
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intensity  is  commonly  5dB  to  lOdB  above  the  active  intensity.  Thus,  care  must  be 
taken  to  insure  phase  matching  throughout  the  real  time  analyzer  if  the  analyzer 
is  going  to  be  used  for  nearfield  measurements. 

The  next  series  of  measurement  investigated  diffraction  which  can  only  be 
identified  with  nearfield  measurements.  Diffraction  was  localized  in  the  nearfield 
of  a  loudspeaker  using  the  reactive  intensity  and  the  imaginary  part  of  the  normal 
impedance.  Diffraction  was  also  studied  at  the  edge  of  a  circular  disc.  A  rigid 
disc  and  an  absorptive  disc  were  contrasted  along  with  two  source  positions.  The 
combined  analysis  with  the  loudspeaker  showed  that  the  diffraction  depended  on 
the  incident  angle  of  the  energy  at  the  impedance  discontinuity  which  caused  the 
diffraction.  More  importantly,  these  measurements  demonstrate  that  to  localize 
diffraction,  the  measurements  must  be  made  acoustically  close  to  a  structure. 

The  final  example  further  shows  the  need  to  use  more  than  one  quantity 
to  entirely  characterize  the  sound  field.  A  plate  like  structure  which  was  being 
excited  by  an  incident  sound  field  was  analyzed.  The  reactive  intensity  indicated 
that  the  structure  was  vibrating  in  a  mode,  but  the  active  intensity  showed  that 
no  power  was  being  radiated.  However  the  measurement  frequencies  were  far 
below  the  coincidence  frequency  of  the  plate,  so  the  vibration  radiates  very  little 
power.  The  reactive  intensity  identified  the  plate  vibration  because  the  energy 
was  confined  to  the  nearfield  of  the  plate. 

All  the  measurements  demonstrated  that  the  analysis  of  intensity  measure¬ 
ments  must  consider  all  available  quantities  along  with  a  detailed  knowledge  of 
the  source  being  investigated.  Table  8.1  summarizes  how  the  active  intensity, 
reactive  intensity,  energy  densities,  and  specific  acoustic  impedance  can  be  used 
to  locate  and  distinguish  between  vibration  sources,  resonance,  and  diffraction. 
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Table  8.1  A  summary  of  the  distinguishing  characteristics  of 
active  intensity,  reactive  intensity,  potential  energy, 
kinetic  energy,  and  real  and  imaginary  parts  of  the 
specific  acoustic  impedance  near  vibration  sources, 
resonators,  and  diffraction. 


Measured 

Quantity 

Vibration 

Sources 

Resonators 

Diffraction 

I 

local  maximum 
points  in  or  out 

local  maximum 
points  in 

not  greatly 
altered 

Q 

local  maximum 
points  out 

local  maximum 
points  in 

points  out 

V 

local  maximum 

local  minimum 

local  maximum 

T 

local  maximum 

local  maximum 

local  maximum 

R'{£«) 

local  minimum 

local  minimum 

sharp  change 

local  minimum 

deep  minimum 

sharp  change 

Most  sound  sources  are  more  complicated  than  the  examples  used  in  this 
thesis.  However,  the  conclusions  can  be  used  to  suggest  a  procedure  to  investigate 
more  complicated  sources. 

Using  the  above  conclusions,  let  us  consider  a  measurement  which  would 
be  made  on  a  motor.  The  motor  and  its  surrounding  structure  have  areas  of 
vibration,  resonance  in  the  motor  enclosure,  and  diffraction  by  the  surrounding 
structure.  These  features  all  contribute  to  the  sound  radiated  to  the  farfield, 
where  the  separate  contributions  of  these  features  cannot  be  identified.  However, 
in  the  nearfield  the  active  and  reactive  intensity  will  locate  each  feature.  Then 
the  features  can  be  classified  as  vibration  sources,  resonators,  diffraction,  or 
interference  effects.  With  such  detailed  information  the  optimum  procedure  to 
reduce  the  sound  radiated  to  the  farfield  can  be  addressed. 
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8.2.  Recommendations  for  Further  Work 

This  thesis  has  been  confined  to  monochromatic  sound  fields  to  establish  a 
clear  understanding  of  the  active  and  reactive  intensity  vectors  for  the  simplest 
case.  However,  the  majority  of  applications  concern  broad  band  and  impulse 
sources.  Also  many  of  the  measurements  in  industry  are  made  in  one-third 
octave  bands  and  octave  bands. 

The  interpretations  of  concepts  such  as  impedance  and  wave  propagation  are 
not  clearly  defined  in  broad  band  sound  fields.  These  concepts  were  studied  in 
this  thesis  by  explicitly  expressing  the  time  dependance  of  a  monochromatic 
sound  field.  Thus,  in  a  monochromatic  sound  field  the  time  dependent 
phenomena  such  as  wave  propagation  and  energy  propagation  could  be  developed 
from  the  time  independent  quantities  which  are  measured.  However,  with  broad 
band  sound  this  mathematical  link  between  the  time  independent  and  time 
dependent  descriptions  is  lost,  because  the  time  dependence  cannot  be  recovered 
after  time  averaging  in  a  broad  band  sound  field.  Therefore,  the  interpretations 
of  the  intensity  vectors  must  be  extended  to  include  multi-frequency  signals. 

Even  at  single  frequencies  the  instantaneous  intensity  is  not  a  plausible  means 
to  describe  the  propagation  of  acoustic  energy.  The  time  independent  active 
intensity  represents  the  wave  front  propagation  and  the  time  averaged  power 
flux,  but  it  does  not  describe  the  path  or  speed  of  energy  propagation  to  the 
farfield.  The  group  velocity  is  used  in  multi-frequency  fields  to  quantify  the 
average  speed  of  energy  propagation,  but  it  is  derived  for  the  time  averaged 
farfield.  Consequently,  the  knowledge  of  the  group  velocity  needs  to  be  further 
developed  to  incorporate  the  time  dependent  process  of  energy  propagation  and 
the  energy  propagation  in  the  nearfield. 
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A  very  interesting  application  of  the  nearfield  acoustic  intensity  measure¬ 
ments  is  to  combine  the  analysis  with  the  structural  intensity  measurements.  By 
measuring  the  energy  in  both  the  structure  and  the  fluid,  the  energy  propagation 
within  the  structure  can  be  understood  along  with  the  coupling  to  the  acoustic 
field.  Such  a  combined  analysis  would  be  invaluable  for  large  structures  such  as 
airplanes  and  ship  hulls.  However,  before  the  analysis  can  be  used  on  complicated 
structures,  simple  panels  with  ribs  and  varying  supports  must  be  studied. 

Nearfield  analysis  of  simple  structures  has  been  recently  published  (Williams, 
Dardy,  and  Finke  1985;  Williams,  Dardy,  and  Washburn  1987);  however,  the 
reactive  intensity  and  impedance  was  not  considered.  Although  the  time  averaged 
power  balance  between  the  structure  and  the  sound  field  is  adequately  described 
by  the  active  intensity,  the  coupling  between  the  structure  and  the  sound  will  be 
strongly  influenced  by  the  reactive  intensity.  Consequently,  the  reactive  intensity 
and  impedance  must  be  included  in  the  analysis 

Throughout  the  measurement  examples  the  distance  of  the  measurement 
plane  to  the  surface  being  investigated  had  the  greatest  effect  on  the  impedance 
measurements.  This  was  especially  important  when  the  intent  of  the  measure¬ 
ment  was  to  represent  the  impedance  at  the  surface.  When  the  measurement 
distance  was  greater  than  ^  the  impedance  was  greatly  affected  by  the  variations 
in  the  acoustic  field  above  the  surface.  Consequently,  the  impedance  on  the 
measurement  plane  was  not  a  good  representation  of  the  impedance  at  the 
structure  surface.  A  greater  number  of  measurements  should  be  performed  to 
establish  the  required  measurement  distance  from  a  structure  radiating  sound 
for  the  measured  impedance  to  be  an  accurate  representation  of  the  impedance 
at  the  structure  surface. 

The  reactive  intensity  is  mainly  concentrated  in  the  nearfield  of  an  acoustic 
source,  so  that  the  power  radiated  to  the  farfield  is  quantified  by  the  active 


intensity  alone.  However,  if  an  operator  of  a  machine  is  in  the  acoustic  nearfield 
then  the  reactive  intensity  may  be  significant.  So  this  energy  needs  to  be 
quantified  in  a  convenient  manner.  Rather  thc.n  expressing  the  output  power 
of  a  source,  a  complex  source  strength  combined  with  a  directivity  function  may 
be  more  useful.  This  would  quantify  the  power  radiated  to  the  farfield  and  the 
reactive  intensity  in  the  nearfield. 

The  pressure  phase  distribution  can  be  calculated  from  the  time  independent 
active  intensity  and  potential  energy.  The  pressure  amplitude  and  phase  as 
a  function  of  space  can  be  used  for  acoustic  holography.  Thus  holography 
measurements  can  be  performed  using  a  single  intensity  probe.  Single  pressure 
microphone  holography  can  be  performed,  but  the  relative  pressure  phase  at  each 
measurement  point  must  be  carefully  monitored  with  a  reference  microphone  or 
with  direct  reference  to  the  vibration  of  the  source.  Using  the  intensity  technique 
no  such  monitoring  is  required. 

A  velocity  sensitive  probe  would  greatly  increase  the  usefulness  of  the 
intensity  technique  because  the  signal  processing  would  be  much  easier  and  a 
larger  frequency  range  could  be  measured.  There  is  one  probe  built  by  Norwegian 
Electronics  which  uses  ultrasonic  doppler  shifts  to  measure  the  particle  velocity. 
However  this  probe  is  large.  Recent  developments  by  Josserand  and  Maeri’eld 
(1985)  showed  the  use  of  a  piezoelectric  material,  PVDF,  for  a  velocity  sensitive 
hydrophone,  which  simultaneously  measures  the  pressure  and  velocity.  This 
material  also  works  in  air.  This  is  a  very  promising  material  and  should  be 
explored  in  greater  detail. 

Measuring  the  particle  velocity  directly  with  a  laser  should  also  be  investi¬ 
gated.  The  laser  only  measures  the  particle  velocity,  but  the  pressure  can  be 
calculated  using  the  equation  of  continuity  which  equates  the  pressure  to  the 
spatial  derivative  of  the  particle  velocity.  This  is  the  direct  analog  of  using  the 
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Euler’s  equation  to  calculate  the  particle  velocity  from  the  spatial  derivative  of 
the  pressure.  The  tremendous  advantage  of  a  laser  based  intensity  system  is 
that  there  is  no  probe  which  would  disturb  the  acoustic  field.  Likewise,  the 
measurements  can  be  made  much  closer  to  a  surface  than  the  two  microphone 
probe.  Further,  it  is  possible  to  simultaneously  measure  three  spatial  components 
of  the  particle  velocity. 
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